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ABSTRACT 


Energy conversion between potential and kinetic energy is considered. Section 2 contains the derivations 
which are necessary to compute the energy conversion for a large region on the basis of vertical velocities and relative 
topography as obtained from a two-parameter model presently used by the Joint Numerical Weather Prediction 
Unit. The energy conversion is divided into three parts: (1) energy conversion due to a mean vertical velocity 
over the region, (2) energy conversion in meridional planes, and (3) energy conversion in the zonal planes. 

Section 3 contains a discussion of the results obtained for the months January and April 1959. The energy 
conversion is positive for each day in both months, but the conversion in the meridional planes has a different sign 
in the two months, being positive in January and negative in April. The pattern of the mean meridional circulation 


is discussed and the frictional dissipation estimated. 


Section 4 describes a procedure for an evaluation of the energy conversion for the different wave numbers, and 


discusses the results for the same two months. 
a linear, adiabatic theory. 


The same section contains a comparison with results obtained from 


Section 5 contains a discussion of the modifications to the results in section 4 caused by the diabatic heating of 
the atmosphere It is made plausible that the maximum conversion found for the small wave numbers by an 
adiabatic computation is greatly altered due to the effects of the heating. 


1. INTRODUCTION 


The problem of conversion between potential and kinetic 
energy has been studied by numerous authors from 
Margules’ [9] original attempt to account for the kinetic 
energy of storms to very recent investigations by Phillips 
[13] and Smagorinsky [16] of the general circulation and 
the energetics of the atmosphere. These latter investi- 
gations are concerned with numerical integrations of 
idealized models of the atmosphere, including, however, 
the effects of diabatic heating and friction. A study of 
ehergy conversions based on data has been restricted in 
the$past mainly because these conversions depend on the 


‘During the preparation of this study the author became aware of a similar study 
conducted by Dr, Saltzman and Dr. Fleischer at Massachusetts Institute of Technology 
Their results for the month of February 1959 agree quite well with those reported here’. 
(Personal communication.) 
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field of vertical motion in the atmosphere. The opera- 
tional use of simplified baroclinic models of the atmosphere 
for short-range numerical prediction has made available 
daily fields of the vertical motion at at least one pressure 
level in the atmosphere. The vertical motion is of course 
computed from the model used for the numerical predic- 
tion and does not include the effects of heating and 
friction. It is computed under the assumption that the 
flow in the atmosphere is quasi-geostrophic, frictionless, 
and adiabatic. Consequently, it should be treated with 
great caution. 

The vertical velocities computed from the method men- 
tioned above seem nevertheless to be of approximately 
the correct order of magnitude and with the correct sign 
as is indicated by the use of the quantity in quantitative 
precipitation forecasts, for instance. It seems therefore 
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worthwhile to use the fields of vertical velocities to obtain 
estimates of energy conversions. An investigation of 
this kind was made by White and Saltzman [23] using 
the vertical velocities from the thermotropic model which 
was integrated by Gates et al. [5] for the entire month of 
January 1953. However, the region covered only ap- 
proximately North America. A later investigation by 
Palmén [10] used a somewhat larger region, but this 
computation was made for only a few days in January 
1956. The Joint Numerical Weather Prediction (JNWP) 
Unit has fields of vertical velocities available twice daily 
since November 1958. The octagonal region used in the 
numerical prediction covers the Northern Hemisphere 
down to approximately 13° N. It is these vertical 
velocities which have been used in the investigation 
reported in the following sections. 

The first problem considered in this paper is the total 
energy conversion within the octagonal region. In con- 
nection with energy conversions it is always interesting 
to divide the total conversion into at least two parts, 
the first being due to mean meridional circulations, the 
other to circulations in the zonal planes. The computa- 
tions were arranged to give this information together 
with the mean meridional circulation. 

A second problem has been to study the energy con- 
version as it appears on the different scales of the atmos- 
pheric motions. 


2. CONVERSION OF POTENTIAL TO KINETIC ENERGY 


The computations which have been made assume that 
the atmosphere is in hydrostatic equilibrium. We are 
further going to consider only the kinetic energy of the 
horizontal flow, because the vertical velocity is small 
compared to the horizontal wind. An equation for the 
change of kinetic energy may then be obtained from the 
horizontal equations of motion in the following way. Let 
us first write the equations of motion with pressure as the 
vertical coordinate: 


ou Oe 

ov Op 

= fu. (2.2) 


In these equations V= V (u,v) is the horizontal velocity; 
w=dp/dt, the “vertical velocity’; ¢=gz, the geopotential ; 
g is the acceleration of gravity; zis the height of an isobaric 
surface; and f is the Coriolis parameter. 

Multiplying (2.1) by wu and (2.2) by v, adding the two 
equations, and defining k= (u?+v’), we obtain 


Ok 


2. 
op V:Vo (2.3) 


ok 
at Vk-+w 


Integrating (2.3) over the complete volume of the atmos- 
phere and defining the total kinetic energy by: 
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where p is the density and S the area over which we inte- 
grate, we obtain 


at ta do v,dldp+ = dldp dSdp 
(2.5) 


In deriving (2.5) we have used the boundary conditions 
w=0 for p=0 and p=po. We have further transformed 
certain area integrals into line integrals using Gauss’ 


theorem (f means a line integral around the closed 
L 

boundery ZL of the region S, while 2, is the outward 

directed normal velocity ) 


The terms on the left side measure the change in kinetic 
energy of the volume, and the transport of kinetic energy 
and potential energy across the boundery. The term on 
the right hand side measures the conversion of potential 
to kinetic energy within the volume (White and Saltzman 
([23]). We shall, with Phillips’ notation, write 


ase 
{P, fies cp (2.6) 


where {P, K} means the conversion from potential to 
kinetic energy. 

The expression (2.6) was evaluated in the following way. 
The vertical velocities available apply at the 600-mb. level 
according to the model approximations, which further 
include the assumption that w has a parabolic distribution 
with zero points at the 1000- and 200-mb. levels. The 
derivative was evaluated by finite differences as 


Oz 1 
dp h (2.7) 
where A is the thickness field between 850 and 500 mb. 
This approximation for 02/dp was used only because h 
was available for each day. 
When the two approximations are introduced in (2.6) 
we arrive at the following formula 
32 
{P, K}=—=] whdS (2.8) 
21Js 
where w now is the vertical velocity at 600 mb. 
Let us next separate w and A into their mean values and 
deviations from these mean values. We define 


h=h+h’ 
(2.9) 
w=o+w’ 
where the mean values are defined as area means: 
(2.10) 
(as. 


a= 


6) 


8) 


nd 
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]serting the expression (2.9) in (2.8) we obtain 


{P, K}= ohS ods. (2.11) 

The first term in (2.11) represents an energy conversion 
due to a mean vertical velocity over the region. Now, 
if S were the complete surface of the sphere we would 
have because 


V-V dSdp=0, (2.12) 
0 s 


i. e., the divergence integrated over the sphere is zero. 
For a restricted region S, as used in these computations, 
there is no assurance that 2=0. We find in general from 
the data that @ has small positive or negative values. 
The first term in (2.11) was computed separately. 

It is furthermore of interest to divide the last term in 
(2.11) into energy conversion due to the presence of a 
mean meridional circulation. We define therefore the 
mean meridional circulation by the following operator: 


L 
( ) dr (2.13) 


where L is the length of the latitude circle and z the east- 
ward directed coordinate. 


Writing now: 


w’ (2.14) 
we obtain 
> 32 _ 32 __ 32 
—{P, K},+{P, K}2t+{P, K}s. (2.15) 


The term {P, K},, the energy conversion in the mean 
meridional circulation, was also computed separately. 

The values of the vertical velocity and the relative 
topography used in these computations were given at the 
grid points of a quadratic grid on a polar-stereographic 
projection (standard latitude 60° N.). Special care has 
therefore to be taken when we perform the numerical 
integrations of the different terms in 2.15). Let us con- 
sider, as an example, the last integral {P, K},. We have 


4 32 
Where (AS), is the area on the earth corresponding to an 


elementary grid square on the map. Performing the 
computations on the map leads then to the expression 


Where m is the map scale factor, (1+sin 60°)/(1+sin ¢), 
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and (AS), is the square of the grid interval at the 
standard latitude, in our case (3.8110°)’m.? Similar 


caution was taken in forming the mean values ( ) end 


( ) and in evaluating the other integrals. 

The three terms {P, K},, i=1, 2, 3 im (2.15) may be 
thought of as conversions due to (1) a mean vertical lifting 
or sinking of the air mass within the region S, (2) circula- 
tions in the meridional plane, and (3) circulations in the 
zonal plane. The two lest terms measure also according 
to Lorenz [7] transformations between zonal available 
potential energy and zonal kinetic energy and transforma- 
tion between eddy available potential energy and eddy 
kinetic energy, respectively. 


3. RESULTS OF THE TOTAL ENERGY CONVERSIONS 


The procedure outlined in the preceding section was 
used to compute the three terms in (2.15) for each day in 
the months January and April 1959. Only the data 
applying at 0000 GMT were used in the computations. 

Let us first summarize the results obtained for January 
1959. Table 1 contains in the first row the total energy 
conversion in the average for the month in the unit kj. 
sec.-'. In the second row the same results are rewritten 
in the unit kj. sec.-'m.~? by dividing by the total area 
of the octagonal region (1.97 * 10"m.*). For comparison 
with other estimates the numbers in the second row are 
converted to the c.g.s. units, erg, gm.~', sec. ~' and written 
in the last row. 

Now, in considering the numbers presented in table 1 
we notice a number of interesting things. The first 
column shows (see (2.15)) that the initial vertical veloci- 
ties indicate a mean downward motion. The mean 
vertical motion was in fact downward on 70 percent of 
the days going into the computation. We shall not 
include {P, K}, in our measure of the total energy con- 
version since this term would be zero if we considered the 
whole earth. It would, incidentally, be an advantage to 
impose the constraint in the numerical forecasts 
with two-parameter baroclinic models over an elmost 
hemispheric region in order to avoid fictitious changes in 
total circulation, es shown by the author (Wiin-Nielsen 
{[21]). If G0, as in these computations, we get a sig- 
nificant energy conversion solely due to the mean lifting 
of sinking, amounting to about 10 percent of the sum of 
{P, K}, and {P, K}s. 

The next important quantity to consider is the energy 
conversion due to the mean meridional circulation. Our 
region is so large that it includes the whole polar region, 
the mid-latitudes, and at least part of the Tropics. The 


TaBLe |.—Total energy conversion computed for January 1959. 


Unit (P, Kh Kh, (P, 
kj. sec.-! —3.14 10" +2.04 K 10" +-28.82 X 10'° 
kj. see.'m.~? —1.59 X 10-¢ X +14.60 X 10> 
erg gm.~' sec.~' —0.159 +0108 +1.46 
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Figure 1.—Vertical velocity in units mm. sec.~! averaged along 
latitudes and in time, for January 1959. The solid curve was 
computed from observed data, the dashed curve from 36-hour 
forecasts. 


adiabatic method of computing the vertical velocity 
makes use in our case of the geostrophic wind. This 
method would certainly break down at and near the 
equator. As the mean latitude of the southern boundary 
of the octagon is about 13° N., we can probably believe 
in the vertical velocities even in the southern portions of 
the grid. The computation for January shows that we 
have a positive conversion of potential to kinetic energy 
in the meridional planes, although it amounts to only 
about 7 percent of the energy conversion in the zonal 
planes. It is very interesting first to see what kind of 
meridional circulation we have in the average for the 
month. This mean meridional circulation was found by 
averaging the vertical velocity for each day in the month 
and then computing the average for the month. The 
result is given in figure 1, where the horizontal coordinate 
is latitude from the pole down to 20° N., while the vertical 
coordinate is the mean vertical velocity in the unit mm. 
sec.-'. Figure 1 shows that the mean meridional circu- 
lation for this month definitely consists of the classical 
three cells with rising motion to the south, the maximum 
downward motion around 35° N., the maximum upward 
motion between 55° and 60° N., and downward motion 
over the polar region. 

The dashed curve on figure 1 indicates the mean meri- 
dional motion for January 1959 for the 36-hour forecasts. 
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The three-cell pattern is still very well defined in the 36- 
hour forecasts, but there seems to be a tendency to pro- 
duce great upward motion in the southernmost portions 
of the region, and perhaps to shift the maximum downward 
motion somewhat to the north by about 5° of latitude. 
Now, with a well-defined indirect cell in the middle 
latitudes, this cell alone would give a conversion of kinetie 
energy to potential energy as was also found by White 
and Saltzman [23] in their computation which uses data 
for January 1953 between 35° N. and 60° N. The net 
positive value obtained in this computation shows that 
the conversion from potential to kinetic energy in the 
two direct cells to the north and the south is sufficient 
to give a positive corversion for the whole meridional 
circulation. It should also be stressed that the main 
part (93 percent) of the conversion takes place in the 
zonal planes. 

The presence of the three-cell meridional circulations 
in the initial vertical velocities is especially interesting 
because these vertical motions are computed from a 
model which does not include the effect of friction and 
diabatic heating. It has been pointed out by Phillips 
[12] in a theoretical study that the three-cell meridional 
circulation is a result of the presence of baroclinic un- 
stable waves. Phillips also computed the intensity of 
the vertical velocity in the meridional cross section. 
In order to make this computation he has estimated an 
average warming of the northern latitudes of about 
0.5° C. day~ in order to balance the yearly average cooling 
due to radiation over the northern half of the Northern 
Hemisphere. In this way he arrives at a maximum mean 
vertical motion of about 3 mb. day~! which, converted to 
our units, corresponds to about 0.5 mm. sec.-' The 
computation of the mean meridional circulation made 
here is certainly in agreement with Phillips’ computation 
as far as the pattern and order of magnitude are concerned, 
although it seems that the mean meridional circulation is 
somewhat stronger in the particular month. 

The total energy conversion for January 1959 amounts 
to about 31 & 10’ kj. sec.~' when we add the contribution 
from the circulations in the meridional and zonal planes. 
As our region is very large we are probably allowed in 
the first approximation to neglect the contributions from 
the advection of kinetic and potential energy into the 
region; see equation (2.5). If this is the case, we have 
that the change in total kinetic energy is equal to the 
conversion from potential to kinetic energy minus the 
frictional dissipation. Including now the frictional dis- 
sipation in equation (2.5) we obtain 


oK 
Spa K}-D (3.1) 


where D measures the frictional dissipation. In a long- 
term mean we have 0K/d0t=0 and the conversion ‘rom 


potential to kinetic energy must balance the frictional 
dissipation; 1.e., 
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D={P'K} (3.2) 
lf we assume that 0K/d0t=0 for the month of January 
1959, we find that D=31 10" kj. sec.-' or about 1.6 
1078 kj. m.~? sec.~', a value which agrees rather well 
with the value estimated by Pisharoty [11], who obtained 
2> 10°° kj. m.~* sec.-! The value obtained from the 
present study is about one-third of the estimate made by 
Brunt [1]: 5 & 107% kj. m.~* see.-' White and Saltzman 
{23| obtained tor the middle latitudes a value closer to 
Brunt’s, while Palmén’s [10] recent computation agrees 
well with the value obtained in this study. 

The integrals in (2.5) may be thought of as sums of 
contributions from different sub-regions of S, although 
the contribution from a subregion is not related directly 
to the change of energy within the region because of the 
contribution from the boundary integrals; i.e., transport 
of potential and kinetic energy into the region. In the 
computation of the total energy conversions for the region 
S, the contribution from each 5-degree latitude ring was 
computed. We could look upon these values as changes 
of the kinetic energy, if there were no transport across the 
boundaries, or simply as measures of the correlation be- 
tween the vertical motion and the relative topography 
(i.e., the mean temperature). The contributions from the 
5-degree latitude bands are plotted in figure 2, where the 
horizontal coordinate is latitude and the vertical co- 
ordinate has the dimension of an energy conversion per 
unit area and unit time. As shown by the figure the 
greatest positive correlation is in the band between 40° and 
45° N. and in the polar regions. In these regions there- 
fore in the average warm air is rising and cold air sinking, 
while this to a much lesser extent is true in the latitude 
bands 55°-65° N. and 20°-35° N. The northernmost region 
of these two represents for this month a region where we 
have a small negative correlation, meaning that in the 
average the warm air is sinking and the cold air rising. 
The curve agrees fairly well with the one given by White 
and Saltzman [23] for the latitude band 35°-60° N. 

The computations reported so far give energy conver- 
sions as computed from initial data. There are therefore 
no forecasts involved except for the very first time step 
which in the present scheme is necessary in order to solve 
for the vertical velocity in the adiabatic equation. This 
procedure could actually be avoided by solving the so- 
called w-equation (see later). From the point of view of 
numerical prediction it is of interest to investigate the 
conversion of potential to kinetic energy as the forecast 
progresses in time. This investigation may be made along 
the same lines as outlined in section 2 replacing w and h 
by the corresponding forecast values. 

The models so far used in short-range prediction have 
been adiabatic and frictionless. For such a model the 
energy equation states that the sum of kinetic and po- 
tential energy is constant. As the models also have been 
hydrostatic, the internal energy per unit column is pro- 
Portional to the potential energy, and the sum of the two 
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Figure 2.—Correlation between vertical velocity and mean tem- 
perature expressed as an energy conversion per unit area and unit 
time as a function of latitude for January 1959. The solid curve 
was computed from observed data, the dashed curve from 36-hour 
forecasts. 


will be called, as is now customary, the potential energy. 
As shown by Lorenz [8] the sum of kinetic energy and the 
available potential energy A is also constant for an 
adiabatic model, and he estimates that the ratio between 
kinetic energy and available potential energy is of the 
order of magnitude of 107'. One important question is 
whether the rate of increase of kinetic energy in a model 
used for numerical prediction is about the same as in the 
real atmosphere. We shall be able to answer this question 
by computing the energy conversion at, for instance, 12, 
24, and 36 hours under the assumption that the region is 
so large that it can be considered closed. For simple 
sinusoidal waves it can be shown that the energy conver- 
sion in a quasi-geostrophic model depends upon the phase- 
lag between the temperature (thickness) field and the 
pressure field. If the temperature field lags behind the 
pressure field, as it usually does in the atmosphere, a 
positive conversion wil! take place, while the opposite is 
true in the reverse situation. Now, it has been noted by 
Thompson [18] that one of the errors in the quasi-geo- 
strophic models is that the phase difference between the 
temperature field and the pressure field decreases too 
rapidly as compared to the real atmosphere. We may 
therefore expect that the energy conversion takes place 
too rapidly in the models, if the observation above is 
quite general. 
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Figure 3.—Vertieal velocity, April 1959. For explanation see 
figure 1. 


It is interesting to notice that this is indeed the case, 
The 36-hour forecast for the month of January has an 
energy conversion which in total amounts to 24 10! kj. 
sec.-' as compared to 3110" kj. sec.~! for the initial 
data, or roughly % of what it should be. An inspection of 
the data for the individual days indicates further that 
there were no days when the conversion computed from 
the 36-hour forecast was greater than the one in the 
corresponding initial data. It may therefore be con- 
cluded that the conversion of potential energy to kinetic 
energy in the quasi-geostrophic model gradually decreases 
as the forecast progresses in time. This fact is also shown 
clearly in figure 2, where the dashed curve gives the 
correlation for the month of January between vertical 
velocity and relative topography in the 36-hour forecasts. 
Although the general characteristics of the curve are the 
same, the average level is lower. 

In order to investigate whether the same general charac- 
teristics would hold for another month the computations 
were carried out for April 1959. They are summarized in 
table 2 and in figures 3 and 4. 


TABLE 2.— Total energy conversion computed for April 1959. 


Unit {P, Kh} {P, K}2 {P, K}s 
kj. sec. —7.72X 10" —2.12X 10!" +21.66 1010 
kj. see.~'m.~? —3.92X 10-4 — 1.08 10-4 +11.00% 10-4 
erg gm.~'sec,~! — 0.392 —0.108 +1.100 


Table 2 shows that we again have a negative contri- 
bution from the term {P, A}, indicating a net downward 
motion over the entire region in the majority of the cases. 
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Figure 4.—Correlation between vertical velocity and mean 
temperature, April 1959. For explanation see figure 2. 


In April it turns out that an average downward motion 
exists on all days except one. 

Contrary to the situation in January we find now that 
the mean meridional circulation converts kinetic energy 
into potential energy. The mean meridional circulation 
still has the three-cell pattern, although some irregu- 
arities now exist in the very high latitudes. The mean 
vertical motion is weaker in April than in January, but 
the extreme values appear approximately at the same 
latitudes in the two months. The main difference 
between the two months is probably that the direct cells 
in the low and high latitudes are less developed in April, 
while the indirect cell in the middle latitudes, although 
somewhat weaker, still is well defined and regular. These 
differences account for the different signs of the energy 
conversion in the meridional plane. However, the energy 
conversion {P, K}, in April is still a rather small fraction 
of the conversion in the zonal planes (about 10 percent). 

If we, for April, assume that the change of the total 
kinetic energy is zero (OK/0t=0), we obtain the result 
that the frictional dissipation for this month amounts to 
about 1.0 10-%kj. m.~? sec.~' as compared to 1.610% 
kj. m.~* sec.~! in January. 

The correlation between vertical velocity and relative 
topography, again expressed as an energy conversion per 
unit area and time, is given in figure 4. We find again 
that the greatest positive correlation is around 45° N ., and 


1.2 3 I 
32 d 
1.0 I 
26 
/ * 24 
° 
t 
\ \ \ A ay 
APRIL 1959 : 
‘ 
\ di 
-12 
-16 
LATITUDE di 
we 
Wi 
mi 
stl 
of 
we 
of 
ga 
th 
CO) 
res 
fro 
the 
the 
at 
wh 
sou 


SepTr (BER 1959 


that vigh correlations exist in the high latitudes although 
more irregularly than in January. 

The 36-hour forecast vertical velocity gives the mean 
meridional circulation drawn as the dashed curve in 
figure 3. The main error in the 36-hour forecasts for this 
month is a shift to the north amounting to about 5° of 
latitude for the maxima and minima, and again we note 
that (he forecast vertical velocities in the low latitudes are 
somewhat larger than observed. This difference is proba- 
bly due to erroneous boundary conditions. 

The conversion indicated by the 36-hour forecasts is 
again somewhat smaller in April than the one given by the 
initial data. While the initial data give a conversion of 
about 20 & 10’ kj. sec.-' within the octagonal region, we 
obtain only about 17 X 10" kj. sec.~! for the 36-hour 
forecasts, or 85 percent of the value given by the initial 
data. 


4, ENERGY CONVERSIONS ON DIFFERENT SCALES 


It has been known for a long time that the atmosphere 
is not equally effective in releasing potential energy on 
different scales. According to baroclinic instability theory 
we find certain scales (4000-6000 km.) which are the most 
unstable. It is also known that the unstable baroclinic 
wave converts potential to kinetic energy. The results 
mentioned above have been obtained from theoretical 
studies using a linear perturbation theory. An evaluation 
of the energy conversion on the different scales is therefore 
worth while doing, first of all to get a check on the validity 
of the linear perturbation theories, and secondly to investi- 
gate whether other energy conversions not treated by the 
theory could be important. Saltzman [14] has recently 
considered the energy equation in the wave-number 
regime. We shall here be interested in only the direct 
conversion of potential to kinetic energy. 

Let us return to the equation (2.8) giving the conversion 
from potential to kinetic energy: 


iP, | has, (4.1) 
s 


We shall in the following consider only a sub-region of 
the complete octagonal region. For convenience we divide 
the polar-stereographic map into rings having the center 
at the north pole. The width of the rings is one grid 
increment on the map. With this arrangement we have: 


whdS= | uhR con (4.2) 
0/0 


where PR is the radius of the earth, considered as a sphere, 
dislatitude, Alongitude, and y the measure of length in the 
south-north direction. We may also write (4.2) in the 
form: 


whdS=R ” Ky) (4.3) 
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with 
0 


Let us next write w(y) and A(y) as Fourier-series in the 
form: 


+3 {a,(y) sin (nX)-+b,(y) cos (nd) } 

> (4.5) 
hOy)=Aoly) +24 {A,(y) sin (nd) + B,(y) cos (nad) } 


Inserting (4.5) into (4.4) and using the orthogonality of 
the functions cos (nd) and sin (nd) we obtain: 


I (y) =2mao(y) Ao(y) +7 {an(y)A,(y) +6, (y) B,(y)}. (4.6) 


The formulae for the Fourier coefficients are: 
1 2r 
w(A,y)dr 


2r 
a= w(A,y) sin (nA)ddA, (4.7) 


by) =~ cos | 


and corresponding expressions for Ay(y), A,(y), Baly). 
Let us next consider the evaluation of the complete 
integral. We have 


y) co dy= (4. 


where dS now is the distance on the map. Evaluating the 
last integral by finite differences we obtain 


I(y) cos ¢ dy= &) I, cos o(1+ sin @) (4.9) 
where we have introduced the expression m= (1+sin @)); 
(1+sin¢) and also the symbol j for the counter of the 
rings. With the present JNWP grid jno;=27; the counter 
) is considered to increase from the North Pole toward the 
equator on the map. 

It is now convenient to rewrite the original expression 
in the form: 


N 
{P, K }={P, K} {P, Kb (4.10) 
n=! 
where 
32 AS 
~ 21 (1+sin do) 
> (4.11) 
(1+sin 
j= 


| 
— 
| 
| 
| | 
) 
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The Fourier coefficients were computed using a proce- 
dure developed by G. ’Arnason, formerly at JNWP. 
The arrangement of the grid points in the octagonal 
grid is not directly suited for a Fourier analysis along 
latitude circles. A short description of the procedure is 
therefore necessary. 

Each of the zonal bands with a width AS is divided 
into zones of 10° of longitude (AA=10°). Exceptions 
to this are made for the high latitudes. In each zone k 
we define 


N 
(asin n (a sin nd), 

(4.12) 
(acos n (a cos nd); 


t=1 


N, is the number of grid points in the zone, k. The 
Fourier coefficients are then computed using the formulae 


36 


The increments AX are as follows: 


j>9, AA=10° 

9>j>5, 

5>j>3, AN=30° 

3>j>0, Ax=60° 

where j, as noted earlier, is the counter of grid distances 
from the pole. 

The computations were again made for each day of 
January and April 1959. The energy conversion was 
averaged for the month to get a picture of the mean 
conditions. The results are shown in figures 5 and 6. 

For each day the Fourier analysis was made up to wave 
number 15 to be on the safe side. In the middle latitudes 
(p=45° N.) n=15 would correspond to a wavelength of 
1800 km., which probably is about the smallest wavelength 
we can hope to analyze with any accuracy with our 
present aerological network. It turns out, as can be seen 
on figures 5 and 6, that the energy conversion is negligible 
for n>11. 

The spectra show for both months two rather broad 
maxima for n=2 and n=6 in January and n=2 and n=7 
in April. The maximum corresponding to n=6 or 7 is 
clearly connected with the most unstable baroclinic wave. 
Converting again to wavelength in the middle latitudes 
we find that n=6 and 7 correspond roughly to wave- 
lengths of 4700 and 4000 km., which coincides well with 
the most unstable waves as predicted from the linear 
perturbation theory. 

It is furthermore predicted by linear perturbation 
theory that all waves with a wavelength shorter than 
about 3000 km. for the middle latitudes should be stable. 
This coincides rather well with the abrupt cutoff in the 
spectra for n=10, which corresponds to L=2800 km. 

These results support strongly the results of the linear 
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Figure 5.—Energy conversion as a function of wave number 
averaged in time, for January 1959. The horizontal coordinate 
is number of waves around the hemisphere, while the vertical 
coordinate is energy conversion per unit time north of about 15° N. 


theory. It is in fact not too difficult to see that an adia- 
batic linear theory must give a result of this nature. 

Let us consider a simple sinusoidal wave pattern in the 
atmosphere. The stream function at 600 mb. and the 
thermal stream function for the layer 800 to 600, or 600 
to 400 mb. will be given in the form 


Yo=—Uy+A sin mr 
(4.14) 
y=—U’y+B sin (mr+a) 


where UW’, and UW’ are the zonal winds considered as con- 
stants, A and B the amplitudes, m=2r/L the wave 
number, and a the phase-difference between the thermal 
wave and the stream function. Note that @ is positive 
if the thermal wave is behind the pressure wave. 

We want in the following to find the energy conycrsion 
for the wave pattern. In order to do this we need the 
vertical velocity field. This can be found from the 
w-equation which for the simple baroclinic model under- 
lying the computations here takes the form 


Fi 


oV 


i iL i i i 4 i iL i iL iL iL 


'23 86 6 8 9 0 8 4 16 


Ficure 6.—Energy conversion, April 1959. For explanation see 
figure 5. 


9 2 
(4.15) 


For the details in this derivation the reader is referred to 
a paper by the author [21], in which it also is shown that 
a solution to (4.15) may be written in the form, 


— (Bv’ (4.16) 
m*o 


for the simple waves given by (4.14). In (4.15) and (4.16) 
o=—ad\|n 6/Op is a measure of static stability, P=40 cb., 
and »,=dy./dx, are the meridional velocities 
in the waves given by (4.14). 

The energy conversion per unit area may now easily 
be computed by inserting the expression for w, and y’ in 
the formula 


As there is no variation in the y-direction we may take 
the area S as a unit band; ie, S=1 - L=2 x R cos ¢. 
We obtain then 
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Ficgure 7.—Energy conversion as a function of wave number 
according to linear adiabatic theory keeping all parameters except 
the wave number constant. 


2m 


m 


Expressed in this form we find that the energy conver- 
sion depends upon (1) the zonal thermal wind, (2) the 
amplitude of the meridional wind components, (3) the 
phase-lag between the thermal field and the stream func- 
tion, and (4) the scale (wave number). 

If we next introduce the number of waves n around the 
hemisphere into (4.18); i.e., m=n/R cos , we obtain: 


maz Unar Sin ah? cos 2 ps 
n ‘ot pr ° Re 


(4.19) 


5{P, K}= 


Several investigations, the first one by Charney [2], 
have shown that the meridional components (x, v’) will 
be greatest for the baroclinic waves (n=6 or 7) and for the 
ultra-long waves (n=2 or 3). The wind components 
above in (4.19) will therefore tend to give a maximum 
for the wave numbers mentioned above. Beside this 
there is a definite influence of the scale of the motion as 
expressed by the last factor in (4.19). If we for a moment 
assume that 7 ner and v,,, are the same for all waves 


| | 
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we can investigate the influence of the scale alone. In- 
serting typical values (U’=5 m. sec.~!, 02 maz=Umar= 
10 m. sec.~', ¢=45° N., sin a=107', o=4) we can plot 
the expression (4.19) as a function of n. This is done 
in figure 7, where we also have multiplied by the area to 
obtain a unit comparable with the units in figures 6 and 7. 
We obtain a spectrum with a broad maximum around 
n=8. It is not too difficult to imagine that this spectrum 
would change into one with a structure similar to those in 
figures 6 and 7 when the standard values of v2 mar ANd Upaz 
were replaced by observed values. As mentioned before 
this would tend to give a maximum for n=2 or 3 and n=6 
or 7. We may therefore say that linear adiabatic theory 
accounts for the gross features of the spectra obtained for 
January and April, although there are too many factors 
mar; Umar, a) to make a close comparison. 

Returning now to the nature of the spectra, we con- 
cluded that the maximum appearing for n=6 or 7 could 
be ascribed to the unstable baroclinic waves. In the pre- 
ceding paragraphs we also pointed to the possibility that 
the maximum for n=2 or 3 could appear because of the 
rather great amplitude in the meridional winds for these 
wave numbers, if we made a comparison with an adiabatic, 
linear theory. The following question may now be asked: 
Would the energy conversion spectra be radically different 
if external effects (heating, friction) could be taken into 
account? Suppose for a moment that the spectrum 
would be changed only slightly. This would mean that 
the ultra-long waves (n=1, 2, and 3) would be self-main- 
taining, having an energy conversion which takes place on 
the same scale, which then would balance the frictional 
dissipation. If on the other hand, heating, for instance, 
will produce a system of vertical motions which are cor- 
related with temperatures in such a way that the energy 
conversion is greatly reduced, we will have to postulate 
another mechanism responsible for the maintenance of the 
kinetic energy of these long waves against frictional dissi- 
pation. One possible mechanism would be a transfer of 
kinetic energy from higher to lower wave numbers through 
non-linear interaction, a mechanism which certainly must 
be operating in the atmosphere. It has in fact been shown 
by Fj@rtoft [4] that a transfer of this nature takes place 
in a two-dimensional, non-divergent fluid. The quasi- 
barotropic character of the atmosphere makes it likely 
that the same mechanism to some extent is operating in 
the real atmosphere. 

A definite answer to the question stated above is not 
easy to give, mainly because the distribution of heat 
sources and sinks is not known with any great accuracy. 
The present knowledge should, however, be sufficient to 
discuss at least the order of magnitude of the effect. The 
following section will deal with this question. 


5. ON THE INFLUENCE OF DIABATIC HEATING 
ON THE ENERGY CONVERSION 


In this section we shall first estimate the vertical veloc- 
ities due to a reasonable distribution of heat sources and 
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sinks. Next, we shall use the distribution of the vertieg] 
velocities to estimate the correction in the energy cop- 
version. 

Very little is known about the distribution of the heat 
sources and sinks in the vertical direction. We shall in 
the following consider a component of the heating fune- 
tion dQ/dt, which is the heating per unit mass and unit 
time prescribed by the formula: 


dQ sinmr. (5.1) 


The main assumption made in (5.1) is that the vertical 
variation of dQ/dt can be assumed to be a power function 
of pressure. The exponent 6 determines the decrease of 
the heating with height. Especially on the large scale it 
is reasonable to assume that the main part of the heating 
and cooling is due to interaction between the air and the 
underlying surface. The heating function ought therefore 
to decrease with height. The parameter 6 has in the fol- 
lowing computations been set equal to 3, which gives a 
rather rapid decrease with height, but the value of 6 is 
open for discussion. 

In (5.1) 7 measures the maximum intensity of the heat- 
ing at the surface (p=po). Some recent computations of 
the vertically averaged heating makes it possible to de- 
termine r. The distribution of heating and cooling pre- 
sented by Staff Members, Academia Sinica [17] shows in 
January a well-defined pattern on the ultra-large scale 
(essentially two waves) with an amplitude approximately 
10-* cal. gm.~! sec.~! for the vertically averaged heating. 
This corresponds to a total heating of 0.4185 kj. m7? 
sec.-' The total maximum heating obtained from (5.1) is 


Equating these two we find that 
r~~0.17 kj. t.~! sec. (5.3) 


If, on the other hand, we take the values used by 
Smagorinsky [15] we find that he considers a value of 0.3 
cal. em.~? min.~'=}X0.4185 kj. m.~? sec.~' as a repre- 
sentative value in winter, although an overestimate in 
summer. As seen, this value is only half the value com- 
puted by the Staff Members, Academia Sinica, and would 
correspond to r~0.08 kj. t.~! sec.-' Our computations 
are, however, linear in r, and it will be easy to find the 
corrections due to the assumed intensity of the heating. 

The author, in collaboration with Dr. N. A. Phillips and 
in connection with other problems, has recently made 4 
computation of the heating for the month of January. li 
is hoped that these computations can be described in detail 
later. It suffices here to say that the heating was com- 
puted from the thermal vorticity equation in the st& 
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Figure 8.—The solid curve is the profile of the mean relative topeg- 
raphy along 50° N. for the normal maps of January given as 
departures from the mean value (left scale). The dashed curve 
gives the latitudinal distribution of heat sources and sinks eom- 
puted from the same normal maps (right scale). 


tionary case in a two-parameter model paying attention 
to the influence of mountains and surface friction. The 
two-parameter model that was used was carefully con- 
structed using the scheme described by Eliassen [3]. The 
normal maps prepared by Jacobs [6], Wege [20], and U.S. 
Weather Bureau [19] for the surfaces 1000, 850, 700, 500, 
300, 200, and 100 mb. were used to define the average 
height field z and the thermal field. The heating derived 
this way for 50° N. is reproduced in figure 8 together with 
the mean temperature field h. 

The result of the computation is in many respects similar 
to the one obtained by Staff Members, Academia Sinica, 
although the maximum heating and cooling in our compu- 
tation appears a little more toward the west than in theirs. 
The maxima are of the same order of magnitude. 

The vertical motions produced by the heating and cool- 
ing may be found from the w-equation. Including the 
diabatic heating in the thermodynamic energy equation 
we arrive at an w-equation of the form: 


R11 
(5.4) 


We are here interested in only the vertical motions due 
to heating and cooling and shall consequently consider 
the following equation: 


2 
f2 %w_ Rl» (4 (5.5) 


In the solution of (5.5) we shall assume that o=a/p’. 
For the goodness of this approximation see a paper by the 
author [22]. A solution to (5.5) may now be written in the 
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form 
w(r,p)=F(p) sin mr (5.6) 


where Fp) will have to satisfy the equation 


r 


Solutions to the homogeneous part of this equation are 
of the form p*. Inserting this function in the homoge- 
neous equation corresponding to (5.7) we find that @ has to 
satisfy the equation 


—a— (5.8) 
giving the two solutions 


1 
(14¥1+4am?/f2). (5.9) 


We find further that F(p)= Bp* is a solution to the non- 
homogeneous equation provided 


R 
=$+1, B= — 5.10 
The complete solution may therefore be written: 
F(p) = Bp p+ (5.11) 


The two integration constants C, and C, are determined 
from the boundary conditions w=0 for p=0 and p= pp, 
giving C,=0 and 


C,=— (5.12) 


The complete solution is therefore 


F(p) = (5.13) 


F(p) gives the vertical distribution of the vertical 
velocity, when it is a maximum. This distribution is 
given in figure 9, where we have used the following 
parameters: fo=10~* sec.~', a= 10*m.? see.~*, r=0.08kj.t.~! 
sec.-', (corresponding to two waves 
around the hemisphere). 

Our next problem is to find the energy conversion due to 
this vertical velocity. The most critical feature here is of 
course in which way the temperature (thickness) patterns 
are arranged relative to the heat sources and sinks and 
therefore also to the diabatic vertical motion. Let us 
write the component of the thickness field as 


h=h, sin (mr+ A) (5.14) 


where A, is the amplitude of the thickness field and A 
the phase-lag between the heat source and the thickness 
field. Note that A is positive if the thickness field is 
lagging behind the heat source field. The diabatic vertical 
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Fiaure 9.—Vertieal profile of vertical velocity caused solely by 
heating. 


velocity may be written 
w=w, sin (mr + 7) (5.15) 
where w, now is positive. 


The energy conversion per unit area due to diabatic 
heating may then be computed from the expression 


(5.16) 
21 L Jo 
which reduces to 
16 
{P, K} wh, cos A. (5.17) 


In estimating A it is important to note that the compu- 
tations of heat sources and sinks show that the atmo- 
sphere is heated where it is cold and cooled where it is 
warm. This is seen from figure 8 where the heating is 
approximately 180° out of phase with the thickness pat- 
tern. The Academia Sinica computation shows about 
the same although the phase difference here is smaller. 
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It seems therefore safe to conclude that cos A for the 
month of January is negative for the large-scale motion, 
The amplitude h, was computed as the amplitude of the 
second Fourier component in the thickness in figure §. 

If the heat sources and sinks were situated just below 
the thermal troughs and ridges, we would have cos A= 
—1 and the energy conversion could be estimated to be 
m.) 


{P, K} 2), dia» ~—9X kj. m.~? see." (5.18) 
as compared to 


{P, K} 1.4 10~* kj. m.~? see.-' (5.19) 


Although (5.18) may be an overestimate due to the 
assumption that cos A= —1 and further to the uncertainty 
of the intensity and position of the heat sources and sinks, 
it is safe to conclude that the large energy conversion 
found by the adiabatic computation in winter is not real, 
but that the heating and cooling produce a system of 
vertical motion which is correlated with the mean tem- 
perature in such a way that the energy conversion from 
potential energy to kinetic energy is greatly reduced or 
perhaps even takes the opposite sign. 

If this conclusion is right, it also follows that the very 
long waves must receive the necessary amount of kinetie 
energy to balance the frictional dissipation in other ways 
than through energy conversions from the potential 
energy. One possibility is, as mentioned before, that 
non-linear interaction between shorter and longer waves 
transfers energy to the long waves. This apparently 
means that it is necessary to consider the changes in the 
shorter waves in order to predict changes on the larger 
seale. This conclusion, if true, makes the prediction of 
the large-scale motion extremely difficult. 


Another conclusion of importance for even short-range 
prediction is that any baroclinic numerical prediction 
model, which does not include the effect of diabatie 
heating, apparently converts too much potential energy 
into kinetic energy on the very large scale. For this reason 
elone it seems important to incorporate large-scale heat 
sources and sinks in a baroclinic model. 

It should also be mentioned that diabatic heating may 
modify the energy spectra on the smaller scale, although 
the main effect there would be heating due to condensation 
rather than interaction with the underlying surface. The 
greater part of the condensation in winter is connected 
with the traveling waves (n=6 or 8). Here, however, it is 
most likely that the diabatie heating will work in the 
opposite direction because a release of condensation heat 
can be of significance only if a considerable amount of 
moisture is available, and this seems to be the case only in 
a warm air mass. For this scale we may therefore expect 
that the vertical motions produced by the heating are 


- 
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positively correlated with the thickness. A qualitative 
arg ment of this type indicates that the energy conversion 
on the scale n=6 or 8 may be greater than indicated in 
figures 5 and 6. 


6. GENERAL CONCLUSIONS 


On the basis of computations of energy conversions from 
potential to kinetic energy for individual days it was found 
that each individual day gave a positive conversion of 
potential to kinetic energy. The mean conversions for 
the months of January and April were 15.6 10~* kj. m.~? 
sec.~' and 10 107* kj. m.~? sec.~', respectively. Estimates 
of the frictional dissipation based on these computations 
agree fairly well with values obtained by Pisharoty from 
independent data. It is found that energy conversion in 
the meridional planes is small in magnitude compared 
to the conversion in the zonal planes. The computations 
further show that the weak meridional circulation consists 
of the classical three-cell pattern. A difference between 
January and April is found with respect to energy con- 
versions in the meridional planes. The meridional 
circulation in January has a positive energy conversion, 
while it is negative in April. This seems to be connected 
with a less well-defined direct circulation in the low and 
high latitudes in April. 

It is further found that the energy conversion gets 
smaller and smaller as time progresses in the present 
JNWP forecasts, verifying an earlier observation that the 
phase-lag between the temperature and pressure fields 
becomes small too rapidly in these forecasts. The mean 
meridional circulation is maintained in the forecasts 
although there seems to be a shift of the pattern toward 
the north. 

The energy conversion for different wave numbers is 
computed in section 4. It is found that an adiabatic, 
frictionless computation shows two maxima for the 
wavelengths corresponding to 2 and 6 or 7 waves around 
the hemisphere. 

In section 5 it is shown that the maximum, which 
appears for the very long waves, most likely would 
disappear if heat sources and sinks could be taken into 
account. It seems therefore that the very long waves 
will receive the main part of their kinetic energy through 
a non-linear interaction with the shorter waves. 
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CORRECTION 


Montuiy Weartuer Review, vol. 87, July 1959, p. 281: In column 2, 
line 2, the vertical motion value should be +i mm./sec. 
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ABSTRACT 


A 48-hr. forecast for the entire Northern Hemisphere of a barotropic hydrostatic atmosphere is made with the 
“primitive equations.’’ Overlapping Mercator and stereographic grids are used, together with the finite-difference 
scheme proposed by Eliassen. Initial data corresponded to a Haurwitz-type pattern of wave number 4. The 
initial wind field was nondivergent and the initial geopotential field satisfied the balance equation. The compu- 
tations seem to be stable and well behaved, except for two small temporary irregularities. The amplitude of the 
gravity-inertia waves present in the forecast geopotential field is about 1/30 that of the large-scale field. It can be 
shown that this is due to the neglect, in the initial data, of the quasi-geostrophically conditioned divergence field. 
The computational technique itself therefore does not give any unreal prominence to the ‘‘meteorological noise.”’ 
The computational characteristics and stability criterion of the Eliassen finite-difference system are investigated for 
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a linearized version of the equations. 


1. INTRODUCTION 


The so-called “primitive equations” have not been used 
much in numerical forecasting because of two main diffi- 
culties. First, if the initial wind and pressure fields are 
not known accurately, artificially large gravity waves will 
appear in the forecast [3,8]. Secondly, the computa- 
tional stability criterion for these equations requires a 
time step of at most 10 minutes compared to the 40- to 
60-minute time step allowed in the geostrophic system. 
The development of larger and faster computing machines 
is rapidly eliminating the second difficulty. It also seems 
probable that a gradual improvement of the rawin and 
radiosonde network, combined with special analyses of 
the initial data, may go far toward solving the first diffi- 
culty. A stable and accurate computation scheme is then 
all that will be required to take advantage of the more 
faithful reproduction of atmospheric processes which is 
possible with the primitive equations. (The geostrophic 
system not only fails at short wavelengths [2], but also 
loses its special prognostic value at extremely long wave- 
lengths [1]. In addition, certain important effects such 
as the horizontal variation of static stability cannot be 
incorporated into the geostrophic system [10,11].) 

There are two aspects to the design of a good compu- 
iation scheme for the primitive equations: (a) the finite 
difference equivalents of the partial differential equations 
themselves, and (b) the formulation of lateral boundary 
conditions. Eliassen [6] and Platzmann [15] have dis- 
cusse’ the former and have arrived at a finite-difference 
scheme for the primitive equations which is more efficient 
than the type of finite-differences currently used in geo- 


strophic forecasts. However, boundary conditions for the 
solution of the primitive equations by finite-differences 
also require considerable care in their formulation, as has 
been pointed out, for example, by Smagorinsky [17]. This 
problem is greatly simplified if the lateral boundary of the 
forecast region can be placed on the equator, where suit- 
able symmetry assumptions can be imposed on the fore- 
cast variables. Although the equator as a boundary is 
readily fitted into either a spherical coordinate system or 
into coordinates on a Mercator map, both of these co- 
ordinate systems have singularities at the North Pole. 

In an attempt to avoid this problem, the writer has 
suggested the simultaneous use of a Mercator map in low 
latitudes and a stereographic projection in high latitudes 
[13]. However, the computational stability of this sys- 
tem is then too complicated a question to be examined by 
mathematical analysis. A numerical test of the scheme 
has therefore been made and the results are described in 
this paper. The equations used were those appropriate 
to a homogeneous incompressible atmosphere moving 
hydrostatically. The initial wind and pressure (geo- 
potential) fields were defined mathematically, rather than 
being obtained from a weather map. The computations 
are therefore a test only of this method of computation, 
and do not purport to answer the question of whether 
adequate initial data can be defined for real forecasts. 


2. THE EQUATIONS OF MOTION 


The horizontal equations of motion in spherical coordi- 
nates, assuming hydrostatic balance and _ neglecting 
friction, can be written 


| | 
59 
he 
ith 
her 
en 
us, 


334 
ad cos 6=2a0(\+ 2) sin (1) 
cos 8 OX 
10¢ 
ab=—ad cos 6(A+22) sin (2) 


Here \=longitude, 6=latitude, a=radius of the earth, 
Q=angular velocity of the earth, and ¢=the geopotential 
of an isobaric surface. The dot (°) is the substantial 
derivative: 


(=F +i (3) 


I ere p is the pressure, and it is clear that the independent 
variables which are being used are \, 6, p, and t. Equa- 
tions (1) and (2) do not contain the Coriolis term 204 cos 6 
or the inertia terms 2\@ cos @ and 240. They must be 
neglected for consistency when the hydrostatic approxima- 
tion is used, since their counterparts in the third equation 
of motion have also been neglected. 

Applying these equations to the motion of an incompress- 
ible homogeneous atmosphere with a free ‘surface we find 
that (1) and (2) carry over as written, if ¢ is set equal to 
gz, Where z is the variable depth of the atmosphere. 
(3) becomes simpler by the disappearance of the 0/dp 
operator, since \ and 6 may be taken as independent of the 
vertical coordinate p. The only other equation needed is 
the continuity equation, which for this atmosphere can be 
written 


cos 8) 


Ot 

(Here g, the acceleration of gravity, has been assumed 

constant, and variations in a have also been neglected.) 

We now define the map coordinates for the Mercator 
and stereographic projections as follows: 


X=an, 
Mercator: (5) 
Y=—alnh, 
_ 2=2ah cos 
Stereographic: (6) 
y=2ah sin i, 
h=cos 6(1+sin @)7!. (7) 


Following the procedure outlined in [13], we define 
cos 6=M-X, 


(8) 
M=sec 8, 


so that U/ and V are the horizontal velocity components 
along the X- and Y-axes of the Mercator projection. M 
is the scale factor for this projection. 

For the stereographic projection we define 
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u=—ad cos 6 sin A—ab cos A=m~'z, 
cos 6 cos A— a6 sin A=m~"y, (9) 
m=2(1+sin @)~?, 

so that wu and »v are the horizontal velocity components 

along the z- and y-axes of the stereographic projection, 


m is the scale factor for this projection. (U, V) and (u, ») 
are related by the expressions: 


u=—U sin \— V cos i, 
(10) 
v=U cos \— V sin X. 


The two horizontal equations of motion and the con- 


tinuity equation (4) can now be written in the map co- 
ordinates: 


og oly ol” ML 
e+ vy |+ Vsin of 20 
(11) 


(12) 
ot =—m| = + +05 | + 20 sin @— 
(14) 
(15) 


The Mercator equations (11)—(13) are equivalent to the 
stereographic equations (14)—(16) and also to the spherical 
equations (1)—(4). 

Equations (11)—(13) have a singularity at the pole 
where M=sec 6 becomes infinite. Equations (14)-(16) 
have a singularity at the south pole, where m=2(1+ 
sin @)~' becomes infinite. As described in [13], the finite- 


difference solution of (11)—(16) is to be carried out over 
one hemisphere by applying (11)—(13) on a Mercator grid 
in low latitudes and (14)—(16) on a stereographic grid in 
high latitudes. The Mercator grid should extend from 
the equator to about 43° latitude, with the stereographic 
grid being responsible for the area poleward of this lati- 
tude (see fig. 1). 


The singularities are thereby avoided 
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3. FINITE-DIFFERENCE EQUATIONS 


To solve (11)—(16) we use the finite-difference scheme 
proposed by Eliassen [6] for the primitive equations. To 
describe this scheme in the present context we must define 
four \attices on the Mercator map and four on the stereo- 
graphic. First we introduce the constant space increment 
A to be used on both maps: 


(P=integer). (17) 


The finite difference coordinates will be (p, g) and (7, j): 


X=(p—po)A, p=0, 1, .. ., P+1, 
Mercator: 


Y=(q—q)A, q=0, 1,..., Q+1. 
(18) 


., 


Stereographic: 
.,f+1. 


(j—Jo)A, j=9,1, 


(It will be convenient to take J as an even integer.) The 
Mercator grid is therefore rectangular (P+ 2) x (Q+2) and 
the stereographic grid is square (J+2)x(/+2). The 
numbers (po, go) and (7%, jo) determine the location on the 
map of the origin of each grid. We define the four 
Mercator and the four stereographic grids (denoted by 
A, B, C, and D) by the four sets of values of (po, go) and 
(i, Jo) Shown in table 1. Figure 2 shows the relative 
orientation on the map of the points on the four grids 
which have the same subscripts (p,q) or (7,7). From the 
definition of X in (5) and A ia (17), it can be seen that the 
points p=0 and p=1 on the Mercator grid are identical 
with the points p=P and p=P+1, respectively. These 
extra points are included in the grids merely for con- 
venience in solving the equations. 

The Eliassen computation scheme is obtained by de- 
fining the variables (U, V, @) or (u, v, ¢) on the grids as 
follows. At time t=At, 2At, 34t, . . ., nAt, U (or u) is 
represented at the lattice points of grid B; V (or rv) is 
represented at the lattice points of grid C; and ¢ (or ¢) is 
represented at the lattice points of grid A. At the 
intermediate times, t=1/2At, 3/2At, 5/2At, . . ., (n+1/2)At, 
U (or u) is represented on grid C, V (or v) on grid B, 
and @ (or g) on grid D. Such an arrangement is much 
more efficient than when all quantities (e.g. U, V, and ¢) 


TaBLe 1.—Origin coordinates of the four grids used in the Eliassen 
finite-difference scheme. 


Grid | Po in jo 


1/2+3/4 1/2+1/4 
1/2+3/4 1/2+-3/4 

0 1/2 /2+1/4 
T/2+1/4 


1/2+3/4 


928468—59—~—3 
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Ficure 1.—General orientation of Mercator and stereographic 
grids, shown on a stereographie projection. The stereographic 
grid covers the square dashed area centered at or near the pole. A 
Mercator grid covers the area between the equator (outer circle) 
and the inner circle. 


are defined at each point of a single grid at all time steps. 
(Note the discussion of fig. 5 in section 7.) For con- 
venience in notation we will hereafter use a prime super- 
script (U’, V’, ’, u’, v’, ¢’) to indicate the quantities 
defined at the intermediate times t=1/24f, . . ., (n+-1/2) At. 
For example, @ en will indicate ¢|X=(p—1/2)A, 
Y=(q—1/2)A, t=nAt] (stored on grid A, where po=1/2 
and q=1/2). however, will indicate ¢[X=pA, 
Y=(q—1)A, t=(n+1/2)At] (stored on grid D, where 
Po=0 and g=1.) 

For convenience in writing the finite-difference equiva- 
lents of (11)—(16) we also introduce the notation 6, and 
oc, for the following operators: 


5 pS ” 


Similar definitions hold for ,, 5;, 6), ¢;, and Finally, 
to eliminate unnecessary repetition of letter subscripts in 
the formulae, a quantity such as S,, will be written simply 
as Soo, S, So, ete. 

In the Eliassen scheme, the six equations (11)—(16) 
applied to the Mercator and stereographic projections, 
result in 12 finite-difference equations. Quantities ap- 
pearing on the right side of these equations are understood 
to have the time subscript n in all cases. The symbol / 
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Ficure 2.—Relative positions on a Mercator (or stereographic) map 
of the points on the four grids A, B, C, and D which have the same 
finite-difference coordinates p, q (or 7, j). 


is used for the Coriolis parameter 2@ sin 6. f and the scale 
factors M and m are identified by additional subscripts 
A, B, C, or D denoting the grid on which they are located. 


Voon— Vo n-1=—Moos (%) { dghoo+U 006 


+a Vo0) (6g Vo0) —U At foow [i+ | (20) 
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{i, [ [ Viver 
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(**) 
100 n+ Yoon = 00c A 
{ 17 +o,[(opi1) } —UpAt 


$00 n+ 1 — $00n= — Mina 


The first six of these refer to the Mercator grid, the second 
six to the stereographic grid. [The numbers i and j) have 
the value (//2+1/4) in equations (25) and (29), and the 
value (//2+-3/4) in equations (26) and (28).] 

The order of solving the equations is (19)—(21) and 
(25)-(27), followed by (22)—(24) and (28)-(30). From 
the form of the equations it is clear that they can all be 
solved on the interior points of the grids (p=1,2,3,.. . , 
and @=1,3,3,.. ., 
J). A further inspection will show that the left (<=0 or 
p=0) and top (j=/+1 or g=Q+1) boundaries can be 
computed with equations (21) and (27); the left and bot- 
tom (j=0 or g=0) boundaries can be computed in (19), 
(23), (25), and (29); the bottom and right (¢=/+1 or 
p=P+1) boundaries in (24) and (30); and right and top 
boundaries can be computed in (20), (22), (26), and (28). 
The remaining boundary values (in general, two adjacent 
boundaries for each grid of data) must be computed by 
other means than (19)—(30) in order to regenerate the 
complete grids of data at each time step. The procedure 
that was used bere is described in the next section. 

The system (19)—(30), although written in a form sug- 
gestive of uncentered differences, actually uses centered 
differences. (The truncation error can be expressed as a 
series in A*? and (At)?.) A special starting procedure must 
be used to get the primed variables at time t=1/2At. In 
the test computation described in this paper, an uncen- 
tered step was used to get the initial values of U’’, V’, ete. 
at n=O(t=1/2At). UU’, V’, ete. were initially known at 
‘=0. Their value at t=1/2At(n=0) was obtained from 
(19)-(21) and (25)—(27) by temporarily replacing At by 
1/2At on the right side of those equations, and the second 
term on the left side of those equations by L’”’(t=0), ete. 


4. BOUNDARY CONDITIONS 


We examine first the equatorial boundary condition, 
which is used to specify the variables on the bottom row 
of the Mercator grid (q@=0). If the motion at one instant 
'=t, over the entire sphere satisfies the symmetry condi- 
tions 


(A,B, fo) =A(A,— 8, to), 
A(X, 8, to) =—O(A,—8, to), (31) 
6, to) =(r,— 8, to), 


equations (1)—(4) will preserve this symmetry for all 
future ¢. It is clear that under these conditions computa- 
tions need be made over only one hemisphere, and the 
appropriate boundary conditions at the equator may be 
inferred directly from (31). For the Mercator variables 
U’, V, and ¢, equation (31) may then be written 


U(X, t) =U(X,—Y, t), 
V(X, Y, ) =—V(X,—Y, 0), (32) 
o(X, t) =¢(X,— t). 


As applied to the six Mercator grids of data, we find that 
this implies 
Up (Go=1/2) 


$50 (qo= 1) 
(33) 
Un= Up (qo= 1) 


V p= Vn (Go= 1/2) 
=Pp1 (Go= 1/2) 


(It can be shown from equations (20) that V7}, will always 
be zero.) 

The left and right boundaries of the Mercator grid are 
easily handled by applying the cyclic boundary condition 
that the point (yp=0, g=q) is identical to the point (p= P, 
q=q) and the point (p=P+1, g=q) is identical with the 
point (p=1, g=q). 

On the top boundary (¢g=(Q+1) of the Mercator grids, 
V’, ¢’, and UU’ can be forecast by (20), (21), and (22). 
The variables L’’, V, and ¢ at g=(@+1 cannot be forecast 
by (19), (23), and (24), however. (Note that these latter 
quantities are stored on grids A and C, which have the 
smaller value of gp in table 1, and therefore the most 
northerly position of the four Mercator grids, as shown 
by fig. 2.) In addition, boundary values of the six stereo- 
graphic variables wu’, v’, ¢’, u, v, and ¢ must be obtained 
on the boundaries listed in the third column of table 2. 
(These are the only points on those grids which cannot be 
forecast from the finite difference equations (25)—(30).) 
In [13] a method was outlined for obtaining these boundary 
values by interpolation from the associated grid. For 
example, a Mercator boundary value of ¢ is interpolated 
from the corresponding stereographic grid of ¢ values, and 
vice versa. In the case of the velocity components 
(U, V) and (u, v) the relations (10) must of course be used 
to supplement the inierpolation process, as indicated in 
the last column of table 2. 
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TaBLe 2.—Interpolation of boundary values between the Mercator and 
stereographic grids 


Variable Grid Boundaries needed Interpolated from grid— 
q=Q+l1 Ster: wu’ (C) and ov’ (P) 
Mere. (A) q=Q+1 Ster: (A) 
_| Mere. (C) q=Q+1 Ster: u (B) and (C) 
_..| Ster. (C) i=J]+1 and j=/+1 Mere: U’ (C) and (B) 
Wosaasss Ster. (B) i=0 and j=0 Mere: U’ (C) and V’ (B) 

| Ster. (D) i=/+1 and j=0 Mere: (D) 
O... i | Ster. (B) i=0 and j=0 Mere: U (B) and V (C) 
CC) i=J+1 and j=J+1 Mere: U (B) and V (C) 
i=0 and j=/+1 Mere: (A) 


In order to perform the boundary interpolations in a 
‘neat’? manner it is necessary that the individual Mercator 
and stereographic grids overlap one another to a certain 
extent. To make this statement more precise, let us 
define the sub-boundary points of a grid as those points 
located next to a boundary (i.e. where g=Q on the 
Mercator grids or the points with <=1 or J, j=1 or J 
on the stereographic grids). We require that the 
points g=Q+/ on Mercator grids A and C (for which 
Y=(Q+1/2)4) lie north of the most northerly sub- 
boundary point of the stereographic grids. This latter 
_ point is located at a distance r=y2?+y?=(A/4)[(27 
—3)?+1]'” from the pole on the stereographic projection. 
The boundary values of U’, ¢, and V required in table 2 
can then be obtained by interpolation on the stereo- 
graphic grids without reference to the boundary values 
on those grids. Referring to (5), (6), and (7) we find 
that Y and r are related by the equation exp(— Y/a) 
=r/2a. We must therefore have 


exp[—(2Q+ 1)4/P]< (#/4P)[(2I—3)?+ 1]}'”. 


(34) 


Another constraint which should be satisfied is that the 
northernmost boundary point on a stereographic grid 
(at r=(A/4)[(27+1)?+1]'”) be located south of the 
southernmost sub-boundary point on the Mercator grids 
(located at g=Q on grids B and D, with Y=(Q—1)A). 
The boundary values of w’, v’, 9’, u, v, and ¢ required in 
table 2 on stereographic boundaries can then be obtained 
by interpolation on the Mercator grids without reference 
to any Mercator boundary points. This leads to a second 
inequality: 


exp[—(2Q—2)/P]< 


If Q and J are not large enough to satisfy both (34) and 
(35), the interpolation process becomes more complicated 
and will undoubtedly lead to mathematical instabilities. 

The values of P, Q, and J used in this test computation 
were 


(35) 


P=80, Q=12, [=22. 


This gives a grid increment A=2zxa/P of 500.4 km. The 
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corresponding distances on the earth varied from 5004 
km. at the pole (m=1) and equator (M=1) to a minimum 
of about 350 km. (m~M~1.4) in middle latitudes. | 
should be pointed out that horizontal space differences 
in the system (19)—(30) are taken over the distance 4, 
rather than 2A as is customary at present in numerical 
weather prediction. 

It might be thought that all boundary values on the 
stereographic grids, including those which can be forecast 
by (19)—(24), could be obtained by interpolation. This 
procedure was in fact tried in a preliminary computation, 
(It is logically easier from the point of view of the machine 
program to interpolate all four boundaries on the stereo- 
graphic grids than it is to do only those listed in table 2.) 
However, the results of this preliminary computation were 
quite unsatisfactory compared to the results described 
in this paper. When interpolation was done on all 
boundaries, the flow patterns tended to move at different 
speeds on the two grids, and discontinuities developed 
near the grid boundaries. 

The details of the interpolation process were as de- 
scribed in [13]. Computation of one interpolated boundary 
value took less machine time than did a computation of 
one of the equations (19)—(30) at one pont. Thus, only 
about 5 percent of the total computation time was spent 
on the boundary computations. 


5. INITIAL DATA 


The initial velocity and geopotential fields for this test 
computation were defined by a flow pattern of the type 
treated by Haurwitz [7]. The initial velocity field y was 
non-divergent, and given by the stream function y: 


y= —a’w sin 6+a? K cos sin 6 cos Rx. (36) 


w, K, R, and a (radius of the earth) are constants. As 
shown by Haurwitz, a flow pattern like this will, in a 
non-divergent barotropic atmosphere, move from west to 
east without change of shape with the angular velocity »: 


(37) 


The equations used in the computations, however, are 
not those for a non-divergent atmosphere, but for one 
with a free surface. (37) will therefore only be satisfied 
approximately. The presence of divergence in the baro- 
tropic atmosphere, will, as is well known, slow up the rate 
of progression of the flow pattern, especially for small 
values of the wave number PF [16]. 

In the non-divergent barotropic atmosphere treated by 
Haurwitz, the pressure field (p/p) associated wit! the 
initial flow pattern (36) can be readily determined by 
integration of the equations of motion (1)—(2), using the 
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anguiar phase velocity v to evaluate dy/dt. We replace 
ple by The distribution of ¢ obtained in this way is 
giver) in the following formulae: 


(6) +a?B(@) cos RA+<a°C(6) cos 2RA, 


A (8) = w) + (R+-1)e? 
+ 


2(Q+w0)K 
(R+1)(R+2) 


BO) = 


C0) 
c=cos 6. 


Both (36) and (38) satisfy the symmetry conditions (31). 
@ in (38) is an arbitrary constant which will determine 
the average height of the free surface in the atmospheric 
model being used here. This in turn will determine the 
speed of propagation of gravity-inertia waves and also the 
order of magnitude of the divergence in the model (div 
t). 

The initial distribution of y and ¢ used for the compu- 
tations was that given by (36) and (38), with the following 
values for the constants: 


w= K=7.848 10-® sec. (~0.19), 
R=4 
¢o=9.8 (8X 10°) m.? sec.~?. 


These values for w and K give rise to large velocities, the 
maximum values of ad cos 6 (=—a-'dy/d0) and aol = 
(a cos @)—'dy/OA] being about 99 and 65 m. sec.-'. Figure 
3shows the distribution of the initial height of the free 
surface z=¢/g. The total variation of 3.5 km. in z is sev- 
eral times as large as the typical variation in the height 
of the 500-mb. surface in winter. 

It is clear from the way in which this initial ¢-field was 
determined, that y and ¢ together satisfy the so-called 
“balance equation” [3]: 


where 
=sin (20 a cos 06)” 
A,=a~ sec =—y-V(ad cos 8), (39) 
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Ficure 3.—Initial distribution of the height of the free surface, 
shown on a stereographic projection. Only one octant is shown, 
the pattern repeating in the other three octants of the hemisphere. 
The outer circle is the equator. Isolines are labeled in km. 


(V is the horizontal gradient operator on the sphere.) 
The advantage of using ¢- and y-fields which satisfy (39) 
is that O(div y)/Ot is initially zero. According to Charney 
[3], this will result in much smaller amplitudes of the 
gravity-inertia waves than would appear if only the 
geostrophic relation were used to relate the initial ¢- and 
v-fields to one another. 

It is of some interest to examine the vorticity field 
corresponding to (36): 


f+ sin (1+ R)(2+-R)K cos* sin cos Rr 
~f{1.1—1.5 cos‘ cos 4d]. (40) 


(We have here introduced the value R=4 and the approxi- 
mate values w=K~0.1Q). Because of the large value 
of K there are four regions of negative absolute vorticity 
in low latitudes. These regions extend poleward to a 
latitude of about 22° where cost @=(1.1)/(1.5). The 
minimum value of f+ ¢ (reached at @~7.25°) is only about 
—0.092, however. 


6. COMPUTATIONAL STABILITY 


The computational stability of the finite-difference 
equations derived in section 3, with the boundary condi- 
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tions discussed in section 4, is difficult to investigate be- 
cause the forecast equations are nonlinear with variable 
coefficients, and the boundary conditions are not simple. 
Considerable information can be obtained, however, by 
examining the computational stability of linearized ver- 
sions of the forecast equations. To further simplify the 
analysis we make the coefficients constant and consider 
only solutions which are periodic in z and y, as on an 
infinite plane. 

As an analogue to (11)—(13) or (14)—(16), we then treat 
the following simple system: 


ou Ou Ou 


Mo dy (43) 


Here f, uo, %, and ® are constants, while u, v, and @ are 
perturbation quantities. (43) could also include a term 
— + rdb/Oy) = —f(urp—vuy). It also is omitted 
for simplicity. For disturbances of the form exp 7 (yt 
+ar+by), (41)-(43) are satisfied by three solutions for 
the frequency yu: 


— (Atty + bry), 


(44) 


u, corresponds to a geostrophic, non-divergent wave, 
while uw. and ys are gravity-inertia waves. 

We now write (41)—(43) in finite-differences, using the 
same basic scheme that was used to get (19)—(24) and 
(25)-(30). We introduce solutions of the form ¥,,= 
exp 7 (ap+q), where a=aA,B=6A, and p and q are the 
finite-difference space coordinates. For convenience in 
the analysis we introduce the following new quantities: 


Upan=U ny BWV yo, 

= Van 40, 


w= (we SiN 5 COS 5 +r cos 5 sin 
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M=(**) sin 


sin 
(45) 


Un, Va, ..., Hi, are now the non-dimensional ampli- 
tudes of the perturbations and are functions of time 
(ndt). We also define the following matrices: 


U, 
H H,, 
(46) 
iF M 
7=—2| -—iF W N 
M 


The Eliassen-grid finite-difference equivalents of (41)- 
(42) can then be written (after some algebra) in the 
following compact form: 


Z. i+ 


GZ, 


An equation for 7 (or 7’) alone is easily obtained from 
these: 


(48) 


Here L is the matrix /+}3G?, where / is the unit matrix. 
L is Hermitian (/;; is the complex conjugate of /;,) and 
therefore has three real eigenvalues \; and three orthog- 
onal eigenvectors ¢;. Further calculation shows that the 
’s are equal to cos &, where £ can assume any one of the 
three values 


sin (3&,)=W, 
sin (3£) = W—[F?+ M2+ N32, (49) 
sin = N22, 


Expanding Z, now as a series in the orthogonal eigenvee- 
tors € 
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we fiid that (48) leads to the following scalar equations 
for the 6;: . 


6; (2 cos bj, +5; =0, j=1,2,3. (50) 


The solutions for 6,, are therefore 


exp n&;), (51) 


where the three values of £; are still given by (49). Com- 
putational stability is achieved by demanding that all é, 
be real. Referring to (49) we see that this requires the 
following inequality to be satisfied: 


(52) 


It should be noted that a weakening of this to permit the 
equality will allow cos 6=—1, whereupon (50) will con- 
tain an unstable solution of the form 6,=n(—1)”. 

Introducing the definitions of W, F, M, and N from 
(45), and taking the worst possible orientation for wp 
and v, the simplified computational stability criterion 
for this Eliassen grid system can finally be written: 


In this formula At and A are the time and space incre- 
ments over which the partial derivatives with respect to 
tand x (or y) are expressed as finite differences. |vo| is 
equal to yU?+°. As is clear from the preceding 
analysis, the satisfaction of (53) will not necessarily insure 
the stability of a computation where the lateral boundary 
conditions are more complicated than the simple ones 
implied by (45). In such cases (53) is best thought of 
as a necessary, but not sufficient, condition for stability. 

Criterion (53) allows a maximum time step of At= 
12.5 min. to be used in forecasting the flow pattern de- 
scribed in section 5. The test computations were made 
with a time step of (1/7) hr.~8.5 min. It took approxi- 
mately 30 sec. on an IBM 704 to compute one time step; 
that is, to solve the 12 equations (19)—(30) at all points 
concerned, and to do the necessary boundary computa- 
tions listed in table 2. A 24-hr. forecast therefore re- 
quired about 84 minutes of computer time. (Checking 
of the results is not included in this figure.) Any further 
increase in computer speeds, say by a factor of 10, will 
certainly make it possible to use the primitive equations 
over an entire hemisphere for even a multi-level baro- 
clinic atmosphere. 

It is clear from (51) and (49) that the linearized finite- 
difference system (47) possesses siz frequencies. The 
conti uous system (41)—(43), on the other hand, possesses 
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only the three frequencies given in (44). It can be shown 
that three of the six frequencies in the finite-difference 
system are similar in form to the three continuous fre- 
quencies in (44), and that the remaining three finite- 
difference frequencies differ only by a reversal in the sign 
of the advection term W. These extra solutions are 
very similar to the “computational wave’’ which is present 
in the conventional way of solving the geostrophic 
vorticity equation [14]. 


7. RESULTS OF THE TEST COMPUTATION 


A 48-hr. forecast was made from the initial wind and 
pressure fields given by (36) and (38). Since this fore- 
cast cannot be compared with either a real atmospheric 
flow pattern or a mathematically known solution, the 
results will be examined only from the following view- 
points: 

a. Smoothness of the fields in space. In particular, 
the agreement between the stereographic and Mercator 
representations in the areas of overlap (see fig. 1). 

b. Smoothness of the fields in time—the question of 
“meteorological noise.”’ 


Figure 4 shows the forecast field of z=¢/g at 48 hours, 
in the area covered by the stereographic A grid. The 
waves have moved about 18° to the east in approximate 
agreement with (37). Of special interest is the agreement 
between the stereographic isolines and the Mercator 
isolines (heavy dashed lines) in the area of overlap of 
the two grids. In general, the two sets of lines are both 
smooth. They agree with one another extremely well 
except for one area near the upper right corner and 
another smaller area near the lower left corner. The 
maximum value of the difference between the two grids 
of z-values in these areas is about 90 meters—about 1/20 
of the maximum difference in z between a trough and 
ridge at the same latitude. The field of z on that portion 
of the Mercator grid not shown in figure 4 was very 
smooth, even in the low latitude regions where /+¢ 
was negative. 

A severe test of the smoothness of the forecast 2-field is 
shown in part A of figure 5. Here the quantity —42;,+ 
+24 at t= 36 hr. is plotted for an 
area centered near the North Pole. (Only the stereo- 
graphic grid covers this region.) There is some tendency 
for a “checkerboard” pattern to appear, but it is not 
very pronounced. 

Smagorinsky [17] and Hinkelmann [9] have made ex- 
perimental forecasts with the primitive equations which 
were not based on the Eliassen type of finite-difference 
grid. In their scheme, only one grid is used (instead of 
the four grids described in section 3), and the geopotential 
and both velocity components are stored at all points of 
this single grid at all time steps. Time and space deriva- 
tives are expressed as centered finite differences over the 
intervals 2 A ¢ and 2 A, much as is done in the usual 
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Ficure 4.—The forecast field of z at 48 hours, drawn on a stereo- 
graphic map, showing only the area covered by the stereographic 
grid A. Isolines are drawn at intervals of 500m. Thin continuous 
isolines are drawn from the stereographie grid-point values (grid 
A). Heavy dashed isolines are drawn from the Mercator grid- 
point values (grid A). The Mercator grid A extended only in to 
the circle. The small dotted areas show the spacing of the 
lattice points on the Mercator and stereographie grids. 
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Ficure 5.—Grid point values of the quantity (—424+ 241 4 + 


2i-1 4+ 21 441 + 2: ¢-1) at 36 hours for an area centered near the pole; 
(A) from the forecast made with the Eliassen grid system, (B) 
from a special forecast made with the finite-difference system used 
by Smagorinsky and Hinkelmann. Units are in tens of meters. 
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way of making numerical weather predictions with, the 
geostrophic model [4]. In order to compare this method 
with the Eliassen method of solving the primitive equa. 
tions, a special 36—hr. forecast was made with the appro- 
priate difference equations from the same initial flow 
pattern. Figure 5 B shows the resulting field of —42,,+ 
Zita t Zinn at 36 hr. from this special 
forecast. The average value of the plotted numbers js 
almost the same, 7.1 in figure 5A and 8.2 in figure 5B, 
but the range in the plotted values is 25 in figure 5A and 
47 in figure 5B. The tendency to a checkerboard pattern 
is very marked in figure 5B. It is clear then that fore- 
casts made this way will be much more irregular than 
those obtained with the Eliassen grid system. 

The main purpose of the computation described in this 
paper was to test the computational stability of the over- 
lapping stereographic-Mercator grids, since this feature 
of the computation was not amenable to the type of com- 
putational stability analysis carried out in section 6, 
Although the results shown in figure 4 certainly indicate 
that the scheme is at least reasonably stable, two small 
temporary ‘“‘wiggles’”’ did appear during the course of the 
forecast. They did not appear until after 24 hours, and 
as shown by figure 4, had practically disappeared again 
by 48 hours. They appeared only near the top boundary 
of the Mercator grid, the stereographic grid point values 
being quite smooth at all times. Figure 6 shows the de- 
tailed structure at 36 hours of the wiggle located near 
4=56°. The other irregularity was very similar and 
located exactly on the other side of the hemisphere, in 
the same part of the wave-like flow pattern. 

Both of the irregularities seemed to be quasi-geostrophic 
in character, with the U’ and V components following the 
geopotential field shown in figure 6. The writer has not 
been able to isolate the cause of these two “errors,” 
which, although small, disfigure what otherwise seens to 
be an excellent computation. Since they appear only on 
the Mercator grids, and in only two of the four waves, it 
is safe to conclude that they do not represent anything 
real, but represent rather some peculiar type of trunca- 
tion error. In this connection it may be important to 
recall from section 6 that the Eliassen grid system does 
contain three false computational frequencies in addition 
to the three physical frequencies. Experience with the 
geostrophic vorticity equation has shown that such false 
frequencies frequently become important near boundaries. 

In order to give an idea of the amount of ‘‘meteoro- 
logical noise” present in the computations, a record of the 
height z at 2-hour intervals is shown in figure 7 for 3 
selected points. (Unfortunately, a record of the forecast 
fields was printed out only every 2 hours=14 time steps.) 
Point I is located near the equator in one of the regions of 
negative f+¢. Point II is located initially in the trough 
at 46° N. near the top of the Mercator grid. (It is one of 
the points in figure 6, where it is marked II.) Point Il 
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is lo ated 175 km. from the North Pole. All three points 
have been purposely chosen in regions of small net changes 
in z ‘except for the last half of curve IT), so that any small 
shori-period oscillations will stand out clearly. Such 
oscillations are indeed present and evidently have an 
amplitude corresponding to a height change of +50 
meters. For comparison, the maximum net 48-hr. change 
in z at any point was equal to +1450 meters-—a value 29 
times as large as the amplitude of the meteorological noise. 
Height tendencies measured over intervals of less than 
about 4 hours would therefore represent primarily ‘“‘noise,”’ 
and not the slower quasi-geostrophic changes. 

This result--the presence of a small but noticeable 
inertia-gravity oscillation--is at first sight contrary to the 
results obtained from “balanced”’ initial data by Charney 
in (3). In Charney’s test computation with the primitive 
equations, no meteorological noise appeared at all when 
the initial wind and pressure fields satisfied the balance 
equation (39). The explanation for this difference is that 
Charney’s initial flow pattern was a stationary wave, while 
the flow pattern used here is not stationary but moves 
slowly to the east. In a barotropic atmosphere a quasi- 
geostrophic wave has, as is well known, a small, but sig- 
nificant divergence field associated with it if the wave is 
not stationary. This divergence associated with the 
geostrophic wave disappears only if the wavelength hap- 
pens to be such that the wave is stationary. Therefore, 
unless the initial wind field also has this small amount of 
divergence, the forecast must contain some high-frequency 
gravity-inertia oscillations (‘‘noise’’?) in addition to the 
low-frequency geostrophic motions. 

From the linearized treatment of the noise problem by 
Hinkelmann in [8], it is possible to estimate the magnitude 
of the noise which is introduced by neglecting in the initial 
data the (small) divergence associated with a moving 
geostrophic wave. If ¢; and c, are the phase velocities of 
the geostrophic and gravity-inertia waves respectively, the 
fictitious gravity-inertia wave will have an amplitude in 
the geopotential ¢ approximately equal to (c,/c2) times the 
amplitude in ¢ of the quasi-geostrophic wave. For the 
example treated in this paper, (¢;/¢2) is about 1/30, giving 
good agreement with the numerically computed ampli- 
tudes in ¢ of the two types of motion. There can be no 
doubt then that the numerically computed noise shown in 
figure 7 is due to the choice of initial data and is not caused 
by the numerical technique. The importance of including 
this geostrophically-conditioned divergence in the initial 
data for the primitive equations has also been demon- 
strated recently by Hinkelmann [9]. 

According to the theory of the geostrophic approxima- 
tion as developed by Monin [12] (contained to some extent 
also in [5]), the second geostrophic approximation to the 
true wind is given by yo, say, where the divergence of y, is 
precisely that divergence which appears (multiplied by /) 
in the usual geostrophic form of _the vorticity equation, 


Ficure 6.— Detailed structure of one of the two temporary irregu- 
larities which appeared on the Mercator grid. The isolines of £ 
(for £=36 hours) are drawn on a Mercator projection, the small 
crosses being the points of Mercator grid A. The top of the figure 
isatqg=Q+1. The heavy dashed isolines represent the Mercator 
analysis and the thin continuous lines an independent analysis of 
the corresponding stereographic grid. Both sets of isolines are 
drawn at intervals of 100 meters, the minimum value isolines in 
the upper right portion of the area being in both cases that for 
8300 m. The grid point marked II is the point corresponding to 
curve IT in figure 7. 


and the vorticity of y, is given by an equation similar to 
the balance equation (39), but with the non-linear terms in 
y evaluated geostrophically. Evidently Charney [3] and 
Hinkelmann [9] have each tested separately the value of 
adding to the geostrophic wind a correction either for the 
vorticity or for the divergence. The elimination of noise 
in both of their results is due to the special choice of initial 
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Ficure 7.—Variation of the height of the free surface at three selected grid points, shown at 2-hour intervals. Units are meters. 


data used by each of them; in Charney’s special test case 
the correction for divergence was unnecessary, as men- 
tioned above, because the wave was stationary, while in 
Hinkelmann’s special case the linear terms in (39) were 
not only much larger than the neglected non-linear terms, 
but the correction for divergence was quite important be- 
cause of the very strong baroclinicity in the zonal flow. 
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Weather Note 


TORNADOES NEAR NAGS HEAD, N.C. JULY 1959 


FRANK B. DINWIDDIE 
Weather Bureau Cooperative Hurricane Reporting Observer, Nags Head, N.C 
(Manuscript received September 15, 1959; revised October 10, 1959] 


On July 10, 1959, at 1300 Est, a tornado was seen to 
form over the mile-wide barrier beach at Nags Head, 
N.C., a little less than a mile from the observer’s house 
(fig. 1). The center of tropical storm Cindy was then 
about 75 miles to the west-northwest, moving northeast- 
ward. Surface conditions at the time were: wind south 
with gusts to 52 m.p.h., sea level pressure 29.79 inches, 
temperature 82° F., dewpoint 72° F., and cloud base level 
about 2250 feet. A loosely formed band of cumulus con- 
gestus was overhead, oriented southwest by south to 
northeast by north, and a heavy rain shower with no 
thunder was just passing over to the north (see fig. 2). 

Suddenly a small cumulus, somewhat rectangular in 
outline, about a mile to the northwest and apparently 
detached from the cloud band, showed unusual though 
not violent agitation as a strong updraft accelerated in 
the body of the cloud and a small horizontal roll developed 
along the nearer edge (figs. 3 and 4). The whole system 
drifted rapidly toward the north-northeast with the gale 
then prevailing. After about 20 seconds the small roll 
cloud elongated southwestward beyond the end of the 
cumulus and a small almost vertical vortex cloud appeared 
near the ground. Immediately thereafter small filaments 
of cloud raced northeastward along the upper edge of the 
roll cloud (fig. 5). This is the only instance I have seen 
or known about of a horizontal roll cloud becoming a 
‘tornado, and I watched it with unbelieving amazement. 
The bottom of this roll cloud rolled toward the parent 
cumulus, which is opposite to the rotation of the typical 
thundersquall roll cloud. 

Several seconds more and the roll cloud had elongated 
into a true funnel, wider than before, inclined upward 
toward the northeast, and joined by a rather sharp elbow 
to the slender nearly-vertical tube at its southwestern end. 
This latter cloud reached to the ground, rotated cycloni- 
cally, and soon developed a cylindrical sheath of flying 
sand around its base, as in figure 6. All these develop- 
ments, from figure 3 to figure 6, were extremely rapid and 
probably did not require more than 30 or 35 seconds. 
Figure 6 is the situation about the time several buildings 
were severely damaged; wind damage was more apparent 
than low-pressure damage. One brick cabin had the roof 
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Figure 1.—Map showing locations of tornadoes and observers. 
Dashed arrows show estimated tracks of tornadoes or damaging 
winds. 

X Location of damage. 

A Location of tornadogenesis 7/10/59 (see also fig. 2). 

B_ Approximate location of tornado 7/25, first photo by Dan 
Morill. 

C Approximate location of tornado 10 min. after B, second 
photo by Dan Morill (fig. 9). 

1 Location of F. B. Dinwiddie observing tornado of 7/10 
(figs. 2-8). 

2 Camera location and tornado bearing 7/25, first photo by 
Dan Morill. 

3 Camera location and tornado bearing 7/25, second photo 
by Dan Morill (fig. 9). 

4 Location and tornado bearing 7/25, Henry Eden. 

5 Location 7/25 of anonymous viewer of tornado. 

6 Location and tornado bearing 7/25 at time of dissolution, 
Mr. Morrison. 

7 Location and tornado bearing 7/25, James W. Pace and 
sketches (figs. 10-17). 
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FicurE 2.—Plan view of position of observer in relation to rain pyguRe 5.—Elongation of roll cloud beyond edge of cumulus and 
shower and cumulus cloud which developed a tornado funnel, appearance of small vertically-oriented vortex cloud near the 
1300 est, July 10, 1959. ground. 


Ficure 3.—Beginning of roll cloud development. 


Figure 6.—The complete funnel. It was at this stage that severe 
damage to several buildings occurred. 
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to Figure 4.—Fully developed roll cloud. Figure 7.—The tornad> moving away toward the north-northeast. 
moved northeastward about 2 feet but left on the walls, nesses in another building described splintering sounds 
“ and part of the southwestern wall blown inward. One passing through the building as if every timber in the roof 
frame cabin was blown partly off its foundation posts and walls was being shattered although no structural dam- 
nd toward the northeast. One ticket booth for an amuse- age was apparent afterward. 


ment center was overturned and one brick chimney was As the tornado passed out to sea, moving north-north- 
blown off. The major damage was the collapsing of the eastward, the funnel became thick at the top, and it was 
front end, walls, and roof of a concrete-block recreation surrounded by several concentric edges of slowly turning 
building. The front faced the west-southwest. Wit- cumulus, looking somewhat like the wrapped leaves at the 
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Ficure 8.—Sketch of the successive thickenings which seemed to move up the funnel just after the stage of figure 6 


Photo by Dan Meril 


Figure 9.—Photograph of funnel about 3:30 p.m. gsv, July 25, 1959, looking north-northwest from Fort Raleigh, Roanoke Island, N.C 
The spout dissipated a minute or two later. 


end of a cigar (fig. 7). It was now joined to the massive 
cumulonimbus formation from which heavy rain was 
falling northeast of the tornado. The tornado was in 
view for about 15 minutes longer. At this stage (fig. 7) 
an important change took place in the action of the tor- 


nado. Small thickenings of the funnel began from the 
surface and grew upward, being cone-shaped with spit 
ridges (fig. 8). When this formation reached about half- 
way up to the cumulonimbus base level, resistance to its 
further ascent developed in the funnel immediately above 
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in th form of bulbous protuberances. It looked as if a 
down (raft from above was opposing the updraft from 
below at this level. A typical sequence of this process is 
show: in figure 8, viewing from left to right. I saw a 
similar situation in a very large tornado-waterspout on 
August 10, 1953, which actually ejected cloud rings from 
protuberances developing about midway up the funnel 
cloud. 

Also on July 10, 1959, tornado-like damage to electrical 
lines and poles was reported from 4 miles southwest of this 
station and about an hour after the disturbance described 
above, but I have not been able to verify any certain 
sightings of afunnel cloud. The next day, July 11, during 
a highly agitated squall cloud condition which produced 
some very strong localized winds 5 miles to the south- 
southwest at Manteo, N.C., one boy reported seeing a 
funne! cloud, but I personally saw none, though I scanned 
the sky with camera in hand. 

Two weeks later, on July 25 at about 1530 xs, a large 
tornado-waterspout was observed near Colington Island, 
about 6 miles west-northwest of this station (fig. 1). 
Although I missed seeing it completely, it was seen by 
many people along the coast and in Manteo on Roanoke 
Island. Figure 9 is the second of two photographs of the 
funnel taken by Mr. Dan Morill from Roanoke Island 
(fig. 1). Very fortunately, the tornado-waterspout was 
sketched by Mr. James W. Pace, my art instructor and a 
man of well-trained powers of observation, as he viewed it 
from near Kitty Hawk (fig. 1). Figures 10 through 17 
are copies of Mr. Pace’s sketches, as he faced about south 
by west from about 5% miles away. 

Surface conditions an hour before the tornado and when 
rain began were: temperature 89° F., dewpoint 75° F., sea 
level pressure 30.07 inches, and wind light air from west- 
southwest. Prior to the onset of rain a band of cumulo- 
nimbus clouds had slowly developed overhead oriented 
west-southwest to east-northeast. During the rain squall 
here, about an hour before the tornado, the wind backed 
to southeast and east, which is unusual, with gusts to 22 
m.p.h. There was no appreciable drop in dewpoint or 
permanent shift of wind after this disturbance, the wind 
afterward averaging about 8 m.p.h., from the southwest. 

Beginning in figure 10 are the shape sequences of the 
tornado seen and sketched by Mr. Pace. Figure 10 shows 
asheath beginning to grow downward around the cloud 
tube, reaching its maximum extension in figure 11. Two 
smal] poorly-formed pendant clouds are on either side of 
the main tornado cloud in figure 12 as it lifted temporarily. 
Rotation in these smaller clouds could not be verified. 
Figure 13 shows what Mr. Pace called the “dumping 
stage,’’ as he thought it was “dumping water’”’ back into 
the Sound. I think it likely that a downflow in the funnel 
cloud was responsible for the downward progress of the 
humps. Figures 14 and 15 show the weakening and termi- 
nal stages of the tornado. Figure 16 shows what Mr. 
Pace thought was a line of very slender but well-formed 
transl cent tubes, which lasted about two minutes and 
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Figure 10.—Tornado cloud as seen by Mr. James W. Pace from 
about 5.5 miles distant near Kitty Hawk, about 1530 est, July 
25, 1959. 


Figure 11.—Maximum extension of sheath around funnel shown 
in its beginning stage in figure 10. 


a pewoawr 
PENDANT 


Figure 12.—Pendulous formations on either side of the funnel 
cloud. 


appeared just after the main funnel disappeared. He 
admitted the possibility that these were thin streaks of 
rain, but did not think they were. The tornado was in 
view about 10 or 12 minutes. 

Figure 17 is Mr. Pace’s impression of a cloud of sand, 
probably whirling, through which he drove several min- 
utes earlier in its formative moments, and which seemed to 
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THE “DUMPING” 
STAGE. 


Figure 13.—A kind of down-flow apparent in the funnel. 


Figure 14.—Weakening of the tornado. 


Figure 15.—Terminal stage of the tornado. 


move very slowly thereafter. This sand cloud was 4 
miles east or northeast of the tornado. No funnel cloud 
was observed near it. 


Figure 16.—Translucent tubes which appeared after main tornado 
funnel had disappeared. 


Figure 17.—Sand cloud 4 miles east or northeast of the tornado 
sketched in figures 10-16. 


Two days later, at 1800 est July 27 and 400 feet from 
where Mr. Pace viewed this disturbance, wind damage 
occurred in a very small area, along Highway 158 at 
Kitty Hawk (fig. 1). The ceiling was blown down in a 
room of a motel, a door blown off, and a boat blown off 
its trailer, while across the road at a service station a 
refreshment stand was blown over a fence toward the 
northeast. The vending machines at the station were 
damaged and its plate glass windows were blown in from 
the southwest injuring three people. Oil cans were 
reported as flying through the air, but whirling motion or 
the sighting of a funnel cloud could not be verified. I 
drove by the spot 30 minutes later, but the rain was 80 
heavy and visibility so poor that I was unaware of the 
damage. Severe lightning was also in progress at the time 
I passed, with the main body of the thunderstorm then 
out at sea to the northeast. 


: Y N \ 
| \ \ 
\ 
ti 
2 
ol 
pl 
he 
pr 
te 
be 
ay 
de 
tic 
eX 
ha 
of 
Aj 
St 
2K 
in 
th 
ey 


BER 1959 


MONTHLY WEATHER REVIEW 351 


MONTHLY PERCENTAGE FREQUENCY OF CEILINGS 2,000 FEET OR 
HIGHER OVER THE UNITED STATES 


M. N. HUNTER, R. I. FOSTER, and H. A. SCOTT 
U.S. Weather Bureau, Washington National Airport, Washington, D.C. 
{Manuscript received June 10, 1959; revised August 28, 1959] 


ABSTRACT 


Maps of the percentage frequency of occurrence of ceilings of 2,000 ft. or higher over the United States, excluding 
Alaska and Hawaii, are presented for each month. A second set of maps shows for each month the change in percent- 


age from the preceding month. 


1. INTRODUCTION 


Low ceilings are frequently the cause of interrupted 
aircraft operations involving either delays or cancella- 
tions of flights. The increase in air traffic and the change 
to jet-type aircraft, which must land at a destination or 
alternate airport with little or no delay, emphasize the 
need for extended weather information for planning pur- 
poses. Government and military personnel, business 
executives, and many airline passengers planning impor- 
tant journeys, need some indication of advance aviation 
weather probabilities to decide upon their mode of trans- 
portation and time of departure. An increasing number 
of requests for information involve aviation weather 
probabilities for periods beyond those for which conven- 
tional aviation weather forecasts can be made at the 
present time. 

If such requests are for periods in the order of 24 to 72 
hours in advance, the meteorologist must consider the 
present synoptic situation and available prognostic ma- 
terial in arriving at some kind of probability guidance. 
Readily available aviation climatological information can 
be used to good advantage as background information in 
such cases. If the request is for an even more extended 
aviation weather probability, the meteorologist must 
depend almost entirely on aviation climatological informa- 
tion. The type of information needed is frequently non- 
existent or inaccessible. 

As an aid to the solution of the problem of providing 
some type of advance aviation weather information, maps 
have been constructed showing the percentage frequency 
of ceilings 2,000 feet or higher over the United States. 
Aireraft operations over the greater part of the United 
States are usually not interrupted when ceilings are above 
2,000 feet, but the problems involving aircraft operations 
inerease as the ceilings fall below 2,000 feet. Therefore 
the 2,000-foot limit was used as the basis for these charts, 
even though it is well-known that terrain, weather, density 


of air traffic, and possibly other things are contributing 
factors in determining the flow of air traffic. 


2. DATA 


The basic data used in construction of these maps were 
obtained from summaries, covering a 5-year period, of the 
hourly observations for 115 airport stations in the United 
States [1]. Supplementary data from [2] and [3] were also 
used with various climatological charts in [4], [5], [6], and 
[7] as supporting material. ‘Topographic maps were also 
used in interpolating over areas where data were sparse. 
Maps, showing the geographical distribution of 2000-ft. 
and higher ceilings were constructed for each month of the 
year (fig. 1) since there is a great variation between 
seasons and also between months of the same season. In 
addition, derivative maps, showing the change from the 
previous month, were constructed for each month of the 
year (fig. 2). 


3. SEASONAL FEATURES 


The following paragraphs describe the most prominent 
seasonal features depicted by figures 1 and 2. 


WINTER MONTHS (DECEMBER, JANUARY, AND FEBRUARY) 


A high frequency of low ceilings prevails over the Great 
Lakes Region and along the western slopes of the Appala- 
chian Mountains from Pennsylvania to Georgia. A 
secondary region of low ceilings occurs from Washington 
State southward over the Cascade Range and Sierra 
Nevada and down through the Rocky Mountains. The 
frequency of low ceilings reaches a peak during the month 
of January over both of these areas with noticeable 
improvement beginning during February. Over southern 
Texas, where a gradual shift in the mean wind to south- 
easterly advects moisture from the Gulf of Mexico, a high 
frequency of low ceiling continues during February. Over 
the remainder of the Plains States and in Florida, rela- 
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Freure 1.—Monthly average percentage frequency of ceilings 2,000 feet or higher. 
tively high, frequency of ceilings above 2,000 feet persists Rocky Mountains may be noted during March with eas 
during the winter season. gradual improvement thereafter. Improvement of the rel 
low ceilings over the Great Lakes Region and the western thi 
slopes of the Appalachian Mountains occurs throughout ine 
Rapid improvement of low ceilings from Washington — the spring months with improvement becoming more rapid alo 


State through the Cascade Range, Sierra Nevada, and during May. Southeasterly advection of moist air ove Ca 
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Figure 1.— Continued. 


eastern ‘Texas and the Southern Plains is evidenced by the 
relatively high frequency of ceilings below 2,000 feet over 
this area, especially during the late spring. A rapid 
increase in the frequency of ceilings below 2,000 feet occurs 
tlong the western sections of the Coastal Range in 
California during April. Elsewhere there is a gradual 


tendency toward a higher frequency of ceilings above 
2,000 feet during the spring. 


SUMMER MONTHS (JUNE, JULY, AND AUGUST) 


Coastal stratus over the western sections of the Coastal 
Range in California prevails throughout the summer. 
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Ficure 2.—Month-to-month change in percentage frequency of ceilings 2,000 feet or higher. 


However, it is interesting to note that the frequency of 
ceilings below 2,000 feet in this area during the summer is 
not as great as that experienced over northern Michigan 
and in part of western Pennsylvania during the winter. 
The remainder of the United States, including Texas and 
the Southern Plains where rapid improvement occurs 


during June, generally has the highest percent of ceilings 
above 2,000 feet during the summer season. 


AUTUMN MONTHS (SEPTEMBER, OCTOBER, AND NOVEMBER) 


Very little change occurs during September except the 
slight increase in frequency of ceilings below 2,000 feel 
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Figure 2.—Continued. 


over northern Michigan, most of Maine, along sections of 
the Appalachian Mountains, and northern Florida through 
‘oastal areas of Georgia. Rapid improvement of the 
‘oastal stratus along California is evident during October. 
The frequency of ceilings below 2,000 feet increases over 


northern Michigan, through the northern Appalachians, 
and from Washington State through the Cascade Range, 
Sierra Nevada, and Rocky Mountains during November, 
with the pace becoming more rapid toward the end of the 
month and through December. 
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THE WEATHER AND CIRCULATION OF SEPTEMBER 1959 


Quasi-Periodic Oscillations in Zonal Index and Centers of Action 


J. F. O'CONNOR 
Extended Forecast Section, U.S. Weather Bureau, Washington,-D.C. 


1. HIGHLIGHTS 


One of the weather highlights of September 1959 was 
the striking of the South Carolina coast by hurricane 
Gracie near month’s end. Due to the timely warnings 
and evacuation of coastal areas, there was no reported 
joss of life from drowning although more than 20 deaths 
occurred in the storm. A few days earlier, September 26, 
one of the most catastrophic typhoons of modern history 
struck Japan near Nagoya. Deaths were reported in 
excess of 4,000, with 1,700 missing, many thousands in- 
jured, and enormous property damage. 

In the United States surprisingly large contrasts in 
weather occurred. Record temperature extremes were 
common. North Platte, Nebr., reported a maximum 
temperature of 102° F. on the 7th, the highest on record 
for so late in the season, and a minimum of 29° F. on the 
10th. At Hartford, Conn., three temperature minima and 
two maxima set new records for individual days. Los 
Angeles had the warmest September of record (average 
temperature 72.8°). The hottest day of the year at 
Buffalo, N.Y., was September 9 when 95° F. was recorded, 
the highest ever reported so late in the season. Bing- 
hamton, N.Y., had the warmest September since records 
began in 1891, while Idaho Falls, Idaho, had the coldest 
September of record. Grand Junction, Colo., had a freeze 
on the 28th, the earliest in 51 years. 

Excessive precipitation occurred this month in many 
areas. For example, Boise, Idaho, reported 2.54 in., 
exceeding all September totals since 1864. Other cities 
reporting record amounts were Kalispell, Mont. with 3.84 
in.; Sexton Summit, Oreg., 2.81 in.; Seattle, Wash., 4.60 
in.; and Oakland, Calif., 3.27 in., the previous record 
being 0.82 in. in 1939. Record snowfall hit areas in the 
lee of the Colorado Rockies. Colorado Springs received 
34.1 in. of heavy wet snow on September 28, 29, and 30, 
with damage to property reported in excess of a million 
dollars. 

By way of contrast this was the driest September on 
record at Lexington, Ky., Milton, Mass., Providence, R.I., 
Oak Ridge, Tenn., and Brownsville and El Paso, Tex. 
Philadelphia, Pa. had its longest consecutive period with- 
out rain, from the 4th to the 27th. 

This was the windiest September on record at many 
Places on the west coast from Sacramento, Calif. north- 
ward to Olympia, Wash. 


2. MONTHLY CIRCULATION AND WEATHER 


The 'argest abnormalities of the Northern Hemisphere’s 
crcula’ on in September 1959 occurred in western Europe 


(fig. 1). A strong anticyclone persisted in the vicinity of 
Great Britain throughout the month, bringing a succession 
of warm, dry, and brilliant days and nights to the British 
Isles, following four consecutive months of unequalled 
dryness, sunshine, and warmth combined. Closely allied 
with the strength of this anticyclone were the intense 
negative height departures from normal to the east, over 
northern Russia, a manifestation of persistent storminess 
and below normal temperatures over eastern Europe 
during the month. 

Over North America, however, there were no unusually 
strong circulation anomalies in the monthly average. 
Figure 1 shows lower than normal heights in the western 
United States and above normal heights in the Northeast. 
These height departures were associated with lower than 
normal temperatures in most of the west and north-central 
sections, while the Northeast, the Mississippi Valley, and 
the Southern Plains averaged warmer than normal (fig. 
2A). Within the month, however, some violent tempera- 
ture fluctuations occurred, and record temperatures at 
both extremes were not uncommon. 

The departures from normal of this month’s precipita- 
tion are shown in figure 2B. September was a wet month 
over most areas except the Northeast and Southwest, in 
line with the prevailing height anomaly pattern which re- 
flected stronger than normal flow of moist air from source 
regions in the Gulf of Mexico, Pacific, and Atlantic. Very 
heavy precipitation occurred in the West north of northern 
Nevada and Utah, with many places in this region report- 
ing this month the wettest or near wettest September on 
record. Most of the heavy precipitation from the Missis- 
sippi Valley westward occurred in the second half of the 
month when the deep mean trough migrated from near 
the west coast to the Continental Divide and became 
stationary there in the latter part of the month. 

Although the average circulation for the month over 
North America was not greatly different from normal, an 
intramonthly upheaval in the broad-scale flow occurred, 
resulting in considerable variability in the circulation and 
weather within the month. This evolution will be con- 
sidered in some detail below. 

3. THE INDEX CYCLE 

An interesting aspect of this month’s weather was the 
oscillation of the speed of the 5-day mean 700-mb. west- 
erlies at temperate latitudes in the Western Hemisphere 
(fig. 3). This was a continuation of a rhythmical varia- 
tion in the strength of the westerlies which started in 
early August with the first substantial weakening of wind 
speeds since June 20. The next low point in the index 
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Ficvre 1.—Mean 700-mb. contours (solid) and height departures from normal (dotted) (both in tens of feet) for September 1959. The 
westerlies averaged near normal, being somewhat depressed in western North America but north of normal in the East. 


(4.8 m.p.s.) was reached at the end of August, after which 
the westerlies increased until the 5-day period ending 
September 7, when the first maximum of September (9.3 
m.p.s.) occurred. This was followed by a decrease to 6.5 
m.p.s. a week later, after which a new high of 10.4 m.p.s. 
was reached on the 26th. Following this peak, another 
decrease set in, which was still in progress at month’s end. 

Figures 4 to 7 are presented to illustrate the sequence 
of 5-day mean circulation and weather anomaly patterns 
associated with this rhythm in the strength of the 
westerlies. 


4. FIRST WEEK—HIGH INDEX 


The high index at the end of the first week was due 
primarily to a deep Low over northwestern Canada (fig. 
4A) which maintained stronger than normal westerlies 
over North America. Five-day mean Lows in September 
have not been observed this far west in Canada in the 
last decade or more, since in this area a ridge normally 
prevails. With the mean Low in this position (see als 
the sea level mean, fig. 4B) the regime was unfavorable 
for anticyclonic activity in Canada. As a result, most of 
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Fievre 2.—(A) Departure of average surface temperature (° F.) 
from normal for September 1959. (B) Percentage of normal 
precipitation for September 1959. (From Weekly Weather and 
Crop Bulletin, National Summary, vol. XLVI, No. 40, October 5, 
1959, and No, 41, October 12, 1959. 


the United States remained warmer than normal except 
for the Far West (fig. 4C). Here cooler than normal 
Pacific air predominated, largely due to stronger than 
normal westerlies along the west coast. The abnormal 
strength of the westerlies can be deduced from the height 
departures in this area (fig. 4A), which totalled 820 feet of 
contour gradient in excess of the normal, between the 
positive center off lower California and the negative 
center in western Canada. 

Figure 4B, which depicts the average sea level flow 
during the first week, shows the predominantly southerly 
flow from the southern Rockies eastward which kept this 
area warmer than normal. Also, the southeasterly drift 
of tropical Atlantic air which favored showery conditions 
tlong the Atlantic Seaboard may be inferred from the 
abnormal location of the Bermuda High near the Mari- 
time Provinces. 

The tracks on figure 4B show that only one High crossed 
the country during the week, and this of Pacific origin. 
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Figure 3.—Time variation of 5-day mean values of the zonal 
index in meters per second, plotted on the last day of the period. 
The index is computed from 35° to 55° N. for the Western Hem- 
isphere. The dashed line represents the variation of normal 
index. Arrows indicate approximate stages of the index cycle 
when amplification, amalgamation, and attenuation of the 
planetary wave train occurred. 


It slowed down over the Northeast and maintained 
easterly flow over the western Atlantic as far north as 
40° N. This flow kept the Southeast wet through much 
of the week (fig. 4D). 

Considerable storminess occurred over western Canada 
with practically no anticyclonic activity. Over the 
United States the only significant cyclonic activity was a 
weak Low which traversed the lower Mississippi and 
Ohio Valleys early in the week, spreading heavy rains 
through much of this area. 


5. SECOND WEEK — AMPLIFICATION AND RETRO- 
GRESSION OF THE LARGE-SCALE WAVES IN THE 
WESTERLIES 


During this week the speed of the westerlies dropped 
sharply in the Western Hemisphere, reaching a minimum 
at week’s end. This was a manifestation of marked 
amplification of the large-scale waves in the westerlies 
over the Pacific and North America. 

Associated with this amplification was an almost com- 
plete phase shift of the major waves over North America 
from the first week. Figure 5A shows the strong upper- 
level High that developed over the Rockies in consonance 
with deepening south of Alaska. This reversal of phase 
might be viewed as a retrogression of the large-scale 
height anomalies over North America and the Atlantic. 
Comparing figures 5A and 4A, it may be seen that the 
Scandinavian positive anomaly retrograded to Iceland, 
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Departure of Average Temperature trom Normal for the Week Ending 
Midnight, 
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Figure 4.—(A) 5-day mean 700-mb. contours (solid) and height ‘ 
a departures from normal (dotted) (both in tens of feet) for Sep- E 
ios tember 3-7, 1959. (B) Mean sea level isobars (solid) (millibars— ' 

hundreds omitted) for September 3-7, with smoothed tracks of 

migratory cyclones (dashed) and anticyclones (dotted) for week 

ending September 7, 1959. (C) Departure of average surface 

temperature from normal (° F.), and (D) total precipitation 

(inches) for the week ending September 6, 1959. (C and D from 

Weekly Weather and Crop Bulletin, National Summary, vol. XLVI, 
, No. 36, Sept. 7, 1959.) 
SEPT. 3-7, 1959 E: 
the positive center over the Maritime Provinces retreated at the same time it deployed cool air into the East as the ing 
to middle North America, and the western Canadian southerly flow of the previous week gave way to northerly flo 
negative center backed off into the northeastern Pacific. flow in this region. This is clearly portrayed by the | 
_§ The amplification of the ridge in western North America mean sea level chart (fig. 5B). The front which ushered en 
occurred almost explosively between September 9 and 10, the cold air into the Eastern States from the 9th to the me 
with simultaneous deepening of the troughs near Kam- 11th dropped temperatures 20°—40° in a few hours, ending me 
chatka and the Gulf of Alaska. During this period the late-season heat and humidity which had prevailed of 
typhoon Patsy (fig. 5B) emerged from the subtropics and for a month or more. This heat wave was climaxed by in 
contributed its energy to the developing Kamchatka record maxima on the 8th and 9th, such as 96° at Sheridan, eX 
vortex. At the same time an upper-level anticyclone Wyo., 103° at Rapid City, S. Dak., 94° at Cleveland, Ay 
from north of the East Siberian Sea drifted southeastward Ohio, and 98° at Rochester, N.Y. Despite the cooling, 10 
and merged temporarily with the intensifying central northern sections still averaged above normal for the th 
Pacific ridge about midweek (fig. 5A) before continuing week due to the large surplus built up in the early partof | 
southeastward to reamalgamate with the amplifying ridge the week (fig. 5C). lig 
in North America. Rains were light over most of the country tis week fir 
The ridge development in the western United States (fig. 5D) since the location of the High in the West pre th 
was accompanied by strong warming over that area while vented moisture from entering most areas. Th: subsid- ex 
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ing motions associated with anticyclonic curvature of the 
low provided an additional desiccating mechanism. 
However, in the Southeast, conditions were quite differ- 
eit. Here cyclonic flow associated with a Low near the 
mouth of the Mississippi River maintained ascending 
notion over the East Gulf States, overrunning the wedge 
of cooler than normal air. This produced excessive rains 
in the East Gulf States and the Carolinas, with totals in 
exeess of 4 inches over a large area. Both Pensacola and 
Appalachicola, Fla., reported over 11 inches, with almost 
\0 inches falling at Pensacola in one 24-hour period on 
the 12th and 13th. 
_Aconiparison of figure 5B with 4B reveals further high- 
lights of the dramatic reversal of the circulation from the 
iirst week to the second week. The cyclone tracks reflect 
the abscnee of cyclonic activity over the United States 
‘xeept in the Gulf of Mexico and in the extreme West. 
Here th» Low from the Gulf of Alaska plunged south- 
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Figure 5.—(A) 5-day mean 700-mb. contours (solid) and height 
departures from normal (dotted) (both in tens of feet) for Sep- 
tember 10-14, 1959. (B) Mean sea level isobars (solid) (milli- 
bars—hundreds omitted) for September 10-14, with smoothed 
tracks of migratory cyclones (dashed) and anticyclones (dotted) 
for week ending September 14, 1959. (C) Departure of average 
surface temperature from normal (° F.) and (D) total precipita- 
tion (inches) for the week ending September 13, 1959. (C and D 
from Weekly Weather and Crop Bulletin, National Summary vol. 
XLVI, No. 37, September 14, 1959.) 


eastward toward the Pacific Northwest as the amplifica- 
tion neared its climax. 

The anticyclone track from the Pacific Northwest to 
the Great Lakes was also a symptom of the extreme ampli- 
fication process going on at upper levels. For one thing 
this High intensified as it moved across the northern 
Plains, and its northward recurvature along Lake Michi- 
gan on the 12th was probably one of the most unusual 
anticyclone tracks on record for this area. The south- 
ward movement of the Gulf of Alaska Low and the north- 
ward movement of the High were both compatible with 
diminishing westerlies at middle latitudes. 


6. THIRD WEEK—AMALGAMATION OF ANTI- 
CYCLONES OVER CANADA 


Early in the third week the large-scale waves reached 
their greatest amplitude. At this time the southeast- 
ward-moving anticyclone from the East Siberian Sea 
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merged over central Canada with the northward-moving 
anticyclone from the western United States. This was 
followed by severing of the anticyclone over central 
Canada as a new branch of westerlies broke across the 
central United States to its south. This evolution pro- 
ceeded in concert with further deepening of the west 
coast trough as the Low from the Gulf of Alaska migrated 
toward the Pacific Northwest. Figure 6A shows the aver- 
age circulation at 700 mb. for period September 15-19, 
with the height departures from normal highlighting the 
blocking character of this stage of the cycle; i.e., a posi- 
tive departure area over central Canada flanked by nega- 
‘tive areas to its southeast and southwest. This period is 
a classic example of the forecaster’s rule that blocking in 
Canada is generally associated with cold and wet condi- 
tions in the United States, shown by figures 6C and 6D. 

The mean sea level chart (fig. 6B) reflects the storminess 
which dominated the Pacific Northwest during this period. 
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Total Precip:tation, Inches, for the Week Ending 


Figure 6.—(A) 5-day mean 700-mb. contours (solid) and height 
departures from normal (dotted) (both in tens of feet) for Sep- 
tember 15-19, 1959. (B) Mean sea level isobars (solid) (milli- 
bars—hundreds omitted) for September 15-19, with smoothed 
tracks of migratory cyclones (dashed) and anticyclones (dotted 
for week ending September 21, 1959. (C) Departure of average 
surface temperature from normal (° F.), and (D) total preeipita- 
tion (inches) for the week ending September 20, 1959. (C and D 
from Weekly Weather and Crop Bulletin, National Summary, Vol. 
XLVI, No. 38, September 21, 1959.) 


The mean Low off the west coast averaged 15 mb. below 
normal, and was largely a reflection of an intense winter- 
type storm which attained a minimum pressure of al 
least 985 mb. near the Oregon coast on the 18th. This 
storm brought high winds and practically all the month's 
record rainfall to northern California. San Francisco te- 
ceived 2.30 inches within 24 hours, more than twice the 
previous September record and 24-hour record. Sal 
Francisco airport also reported high winds of 48 mp. 
on the 18th, 5 m.p.h. higher than the previous record: 

In the East, the week was featured by the southeast- 
ward migration of an unusually intense Arctic antieyelone 
from the Lakes Region across the Middle Atlantic States 
as shown by the track on figure 6B. This track was * 
symptom of the early stages of attenuation o! the extrem 
amplitude of the large-scale flow of mid-month. This 
slowly-moving High produced record-breaking low te 
peratures in the Northeast, from the Great | .akes t New 
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fugland, a marked reversal from the record high tempera- 
lures of the early part of the second week. A few ex- 
amples of the cold accompanying this High are: Lansing, 
Mich., record minima of 32° and 31° on the 16th and 17th, 
spectively; Toledo, Ohio, minimum temperatures in the 
W's on four consecutive days from the 16th to the 19th; 
iid Pittsburgh, Pa., record minima in the 30’s on the 
‘ith, 18th, and 19th. 

Heavy rains this week spread eastward across the 
Rockies to the Great Lakes Region as the southwesterly 
low of moist air that prevailed aloft (fig. 6A) overran 
the cold Arctic air as far east as the Upper Mississippi 
Valley, with Kansas receiving % to 3 inches over most of 
the State. A similar overrunning pattern existed in the 
‘outheast, producing excessive precipitation amounts of 
* 10 inches at seattered places in Florida. 
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Figure 7.—(A) 5-day mean 700-mb. contours (solid) and height 
departures from normal (dotted) (both in tens of feet) for Sep- 
tember 22-26, 1959. (B) Mean sea level isobars (solid) (milli- 
bars—hundreds omitted) for September 22-26, 1959 with 
smoothed tracks of migratory cyclones (dashed) and anticyclones 
(dotted) for week ending September 28, 1959. (C) Departure of 
average surface temperature from normal (° F.), and (D) total 
precipitation (inches) for the week ending September 27, 1959. 
(C and D from Weekly Weather and Crop Bulletin, National 
Summary, vol. XLVI, No. 39, September 28, 1959.) 


7. FOURTH WEEK—ATTENUATION AND PROGRESSION 
OF THE LARGE-SCALE WAVES IN THE WESTERLIES 


The fourth week saw the culmination of this remarkable 
3-week cycle of development and decay of the large-scale 
wave pattern aloft. The attenuation of the major ridge 
over central North America was characterized by an al- 
most complete phase shift over this continent, similar to, 
but in reverse of, the retrogressive shift which accom- 
panied the amplification stage of the cycle in the second 
week. This week’s reversal saw the ridge of the third 
week give way to a trough over central North America 
(fig. 7A). This may be viewed as a progression of the 
Pacific coastal trough eastward in response to increasing 
westerlies. Similarly, the collapse of the block over 
Canada may be viewed as a southeastward progression 
of this feature to the Middle Atlantic States. In terms 
of height anomalies, the negative anomaly off the Pacific 
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Northwest coast moved eastward to the Continental 
Divide, while the Canadian positive center migrated 
southeastward into the eastern United States. 

This evolution set the stage for intensification of cooling 
over the central Rockies, while temperatures in the East 
returned to record warmth (fig. 7C) in many places, due 
to the return of abnormally strong southerly flow. New 
maximum temperature records were set at Akron, Ohio, 
with 88° on the 28th and Pittsburgh, Pa., with 90° on the 
28th and 89° on the 29th, and were equalled at Albany, 
N.Y. on the 24th and 28th with 87° and 83°, respectively. 

The persistent southerly flow during the week from 
eastern Texas tothe Northern Plains and Upper Mississippi 
Valley brought excessive precipitation, high winds, and 
severe local storms to the Midwest, with weekly rainfall 
totals of 1 to 3 inches over a large area (fig. 7D). Local 
precipitation amounts of more than 10 inches in Okla- 
homa, 6 inches in southern Texas and Kansas, and 5 
inches in lowa, Minnesota, Illinois, and Michigan were 
recorded. 

The strength and persistence of the meridionally ori- 
ented polar front associated with the mean trough in the 
center of the country can be assessed by the average 
temperature anomaly gradients for this period in figure 
7C. This strong baroclinic zone bred a number of 
vigorous frontal waves, one of which produced the lowest 
September pressure of record, 989 mb., at Dubuque, 
Iowa on the 26th. This front, overrun by very warm 
and moist tropical Gulf air, spawned several tornadoes in 
Missouri, Wisconsin, Illinois, Kansas, and Iowa, on the 
26th and 27th. Figure 7B shows the mean sea level 
pressure pattern characteristic of the week, with the 
associated cyclone and anticyclone tracks. 

Late in the week the polar front moved eastward as 
cold Arctic air flowed southward into the northern 
Rockies and the high Plains. Overrunning of this wedge 
of cold air produced damaging snowfalls of record intensity 
in eastern Colorado. Pueblo reported a record 14 inches 
of heavy, wet snow which fell on the 29th and 30th with 
heavy damage to trees and shrubs. 

Also near month’s end hurricane Gracie moved inland 
across the South Carolina coast about 11 a.m. Est, 
September 29, a little south of Charleston. Winds were 
estimated near 120 m.p.h. and tides up to 9.7 ft. above 
mean low water in the Charleston area. Heavy rains 
along the storm path (fig. 8) totalled over 6 inches in most 
of South Carolina, 8 inches in western North Carolina, 
10 inches locally in Virginia, 6 inches in western Mary- 
land, 4 inches in West Virginia, and from 1 to 2 inches in 
the New England States. Figure 7D does not reflect 
this rainfall since the hurricane had not moved inland 
until the last two days of the month. 


8. HURRICANES AND TYPHOONS 
ATLANTIC 


Three hurricanes developed in the Atlantic during the 
month. All three developed hurricane winds north of 
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Figure 8.—Selected tropical storm tracks for’ September 195% 
All attained hurricane f°rce except Nora. Dates represent 
approximate 0000 Gar pcsiti~ns. 


20° N., which is unusually far north except for storms 
which originate near the Bahamas as hurricane Gracie 
did. Flora may have originated south of this latitude, 
but Hannah clearly appeared to have originated north 0! 
26°N. All three developed close to the stage of the index 
cycle considered most propitious for hurricanes; i.e., 
times of high index, when the westerlies were farthes! 
north and the subtropical easterlies strongest. Hurr- 
cane Flora developed around the 9th, Gracie on the 22d, 
and Hannah on the 27th, which were times when the 
index was close to its maximum. The development of 
Gracie and Flora appeared to be associated with the well 
known shearing process, in which the southern portions 
of full-latitude major troughs become disassociated from 
the more rapidly moving northern portions near times of 
high index [1]. 

These storms traced widely diverse paths as show? ! 
figure 8. Hurricane Flora followed a very unusual path, 
heading rapidly for the Azores after detection on the th, 
with highest winds reported near 65 knots. ‘lhe eastwar! 
component of this motion throughout its ire durate" 
may be attributed to the sharp drop in the sped and late 
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jude o! the westerlies between September 9 and 12 (fig. 3). 
The lat ter circumstance occurred in consonance with retro- 
gessio!! of the blocking anticyclone from near the United 
Kingdom to a position south of Iceland during this period. 
Asa result the depressed westerlies in the eastern Atlantic 
wept hurricane Flora rapidly northeastward. 

Hurricane Gracie, on the other hand, had its inception 
in the ascendant phase of the index. Thus the westerlies 
were increasing their latitudinal separation from the 
sorm, and it drifted slowly northwestward. Around the 
ime of the index maximum, from the 25th to the 27th, the 
wurricane wallowed indecisively while blocked by a High 
which moved eastward across New England. The hurri- 
cane resumed its northwestward course toward the coast 
on the 27th as the High center moved east of the storm’s 
meridian. This was also about the time hurricane Hannah 
put in its appearance about 1400 miles farther east. Be- 
ween October 1 and 2 both Gracie (now an extratropical 
low), and hurricane Hannah recurved sharply eastward 
as the index dropped below the normal. 

Thus it appears that the Atlantic hurricanes this month 
were spawned when the speed of the temperate westerlies 
and the subtropical easterlies) was above normal, at 
vhich time westward motion prevailed. When the zonal 
index declined below the normal, eastward motion and/or 
reeurvature to eastward occurred. 

PACIFIC 

Inthe North Pacific, 12 tropical cyclones were observed 
this September. In the eastern Pacific, two hurricanes 
leveloped off the west coast of Mexico, and one tropical 
storm occurred east of Hawaii. The first eastern Pacific 
hurricane (fig. 8) was detected near 14° N., 98° W. on 
‘eptember 4, about the time the index rose above the 
wormal. It apparently disintegrated over Baja California 
mn the 11th when the index was dropping sharply. The 


MONTHLY WEATHER REVIEW 


365 


second hurricane (not shown) was detected near 17° N., 
119° W. on the 21st and, after following an irregular path, 
apparently dissipated on the 26th near 20° N., 130° W. 


Tropical storm Wanda, east of Hawaii, was shorilived, 
as was tropical storm Opal and depressions Marge, Ruth, 
and Thelma (not shown). Tropical storm Nora (Sept. 
5-12) moved from east of the Phillipines across Luzon and 
passed near Hong Kong into the Yellow Sea (fig. 8). 
Highest winds were estimated at about 35 knots. 


In addition to typhoon Patsy in the central Pacific, and 
typhoon Louise which developed late in August but hit 
Formosa early in September, two other severe typhoons 
wreaked widespread destruction in the Far East. These 
were typhoon Sarah and typhoon Vera, whose tracks are 
shown in figure 8. 

Typhoon Sarah, first detected on September 11, reached 
typhoon intensity on the 12th, passed over the Ryukyus 
on the 15th, and swept across southern Korea into the Sea 
of Japan. Winds well over 100 knots brought considerable 
destruction to Korea and southern Japan, and heavy 
casualties in the fishing fleets. 

About 10 days later typhoon Vera, a monstrous storm 
in both size and intensity, after originating -near the 
Marianas, moved northwestward and then northward 
along the 135° E. meridian, striking central Honshu on 
the 26th. Winds well over 100 knots plus excessive pre- 
cipitation and tides resulted in enormous casualties and 
destruction, placing it in the category of one of the 
greatest natural catastrophies in the modern history of 
Japan. 
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Weather Note 


A REPORTED SEA LEVEL PRESSURE OF 877 MB. 
C. L. JORDAN 


Florida State University, Tallahassee, Fla. 
(Manuscript received September 8, 1959; revised September 14, 1959] 


Adropsonde observation made in the eye of typhoon Ida 
il approximately 0500 amt, September 24, 1958 showed a 
a level pressure of 877 mb. (25.91 in.) which is believed 
0 be the lowest sea level pressure ever recorded." 2 This 
Sabout 10 mb. lower than the generally accepted mini- 

The value of 873 mb. reported for this typhoon in the 1958 Climatological Data—Na- 
ional Sum mar (p. 93) represents an estimate of the sea level pressure based on the last 
"ported inform:tion from the dropsonde. The pressure-height computation from this 


“unding, whic!) was made 8 hours earlier than the one considered in this note, gave 878 
Mb. a8 the sea level pressure. 


* Spectacular »hotographs of the eye of this storm taken from very high levels have 
Bundgaard, Fletcher, and Smith [1]. 


teen presente: 


mum pressure of 26.185 in. (886.7 mb.) which was observed 
in a typhoon in 1927 by the steamship Sapoeroea east of 
the Philippines, and about 15 mb. lower than the lowest 
Western Hemisphere reading of 26.35 in. (892.3 mb.) 
which was observed at Lower Matecumbe Key during the 
Florida Keys hurricane of September 1935 (cf. [2}). 

The dropsonde observation cited above was made from 
the 700-mb. level by a U.S. Air Force reconnaissance air- 
craft of the 54th Weather Reconnaissance Squadron. At. 
the time of this observation, typhoon Ida was located 
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TaBLeE 1.—Pressure-height, temperature, and dewpoint data for 
standard and significant levels in the dropsonde observation made 
in the eye of Pte, Ida at 0500 amt September 24, 1958. Since 
these data were taken from coded reports, the listed temperature and 
dewpoint values may deviate 0.1° C. from those shown in the original 
records. 


Height (ft.) Pressure Temperat ure| Dewpoint 
(mb.) C.) 
0 877 27.8 26.5 
910 850 26.8 24.8 
835 26.2 23.5 
7 29. 5 15.5 
7 27.8 14.2 
6580 700 29.8 14.2 


near 18.9° N., 135.3° E., or about 600 miles northwest of 
Guam. The temperature, dewpoint, and pressure-height 
data reported in this sounding are shown in table 1. The 
temperature varied over a range of less than 4° C. and the 
maximum value of 29.8° C. (86° F.) was reported at the 
700-mb. level. The 700-mb. surface was reported at 
6580 ft. and the 850-mb. surface at 910 ft. above sea level. 

In the dropsonde observation the 700-mb. height was 
determined by using the pressure and radio altimeters 
aboard the aircraft, and the 850-mb. height and sea level 
pressure were computed using the temperature and 
humidity data measured at flight level by the aircraft 
instruments and at levels below the aircraft by the drop- 
sonde instrument. The sounding under consideration 
was found to be hydrostatically consistent and the vertical 
temperature distribution is realistic, in view of those 
found in other deep typhoons [3]. 

The sea level pressure is relatively insensitive to small 
temperature errors since a relatively thin layer is involved. 
In fact, an error of 5° C. at all levels would change the 
sea level pressure less than 4 mb. in the present case. 
The most likely source of major error in cases of this type 
is in the aircraft altimetry which provides the beginning 
point in the pressure-height computation. Errors of this 
type tend to be systematic for a given aircraft and flight, 
but are normally small on the reconnaissance aircraft 
since calibration checks are made as a routine part of each 
flight. In the case in question, the 700-mb. heights re- 
ported in the vicinity of Guam are consistent with those 
given by radiosonde observations. Any systematic error 
would undoubtedly have been less than 80 ft. which would 
correspond to sea level pressure errors less than 3 mb. 

Perhaps the strongest supporting evidence of the relia- 
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bility of the dropsonde observation under examination js 
the fact that another dropsonde observation made aboy 
8 hours earlier from a different aircraft reported a ge 
level pressure of 878 mb. A third dropsonde mace aboy 
23 hours later reported a sea level value of 886 mb. Sing 
the three dropsonde observations were made from clifferen; 
aircraft, there is little chance that altimetry errors could 
have been important. Even if the improbable assump. 
tions are made that the aircraft pressure-height valyes 
were 150 ft. too low and the reported temperatures were 
3° C. too warm at all levels, the soundings made at 210) 
GMT, September 23 and 0500 emt, September 24 would 
still show sea level pressures lower than the previously 
accepted record value of 886.7 mb. E 

In view of the evidence cited above, it is felt that 
typhoon Ida of 1958 should be accepted as the deepest 
evelone of record. Typhoons with central pressures 4s 
low as those reported in Ida of 1958 are undoubtedly rare, 
However, there have been at least seven other typhoons 
in the past eight years in which the minimum sea level 
pressure was reported to be 900 mb. or lower.’ It is be- 
lieved that typhoon Nina of 1953 would qualify as the 
nearest rival to Ida of 1958 with dropsonde sea level 
pressures of 886 and 883 mb. All of these very deep 
typhoons have contained extremely strong winds and have 
covered very large areas. In the case of typhoon Ida of 
1958, the maximum reported 700-mb. wind speed was 185 
kt. On September 23-24 at the times of the extremely 
low central pressures, typhoon-force winds were reported 
both at the surface and 700-mb. levels at distances as 
great as 275 miles from the storm center. 


ACKNOWLEDGMENT 


The writer is indebted at Lt. Col. Alexander Kouts of 
Air Weather Service for directing his attention to typhoon 
Ida and for providing some of the records used in the 
preparation of this note. 


REFERENCES 


1. R. C. Bundgaard, R. D. Fletcher, and J. R. Smith, ““Typhoon- 
Eye Cloud Patterns as Viewed from Above,” Weatherwise, vol. 
12, No. 2, Apr. 1959, pp. 64-66. 

2. G. E. Dunn, “Tropical Cyclones,” Compendium of Meteorology, 
American Meteorological Society, Boston, 1951, pp. 887-90l. 

3. C. L. Jordan, ‘“‘The Thermal Structure of the Core of Tropical 
Cyclones,” Geophysica, vol. 6, No. 3/4, 1958, pp. 281-297. 


3 These were Marge of 1951, Nina of 1953, Ida of 1954, Virginia, Hester and Lola of 
1957, and Grace of 1958. 


U.S. GOVERNMENT PRINTING OFFICE; 1959 


The 
it irre 
regu 
is “ob 
Ihust | 
itions 
olisiste 
The 
of pol 
forts 
an 
‘igned 
lid no 
ea, 
Wa f 
nes 
bolyno 


=; 
Volume 
| lumber 


